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Abstract 



>< 
^ • The nucleon form factors are calculated using a non-relativistic description in 

terms of constituent quarks. The emphasis is put on the reliability of present numer- 
ical methods used to solve the three-body problem in order to correctly reproduce 
the expected asymptotic behavior of form factors. Nucleon wave functions obtained 
in the hyperspherical formalism or employing Faddeev equations have been consid- 
ered. While a q~^ behavior is expected at high q for a quark-quark force behaving 
like - at short distances, it is found that the hypercentral approximation in the 
hyperspherical formalism [K = 0) leads to a q~'^ behavior. An infinite set of waves 
is required to get the correct behavior. Solutions of the Faddeev equations lead to 
the q~^ behavior. The amplitude of the corresponding term however depends on 
the number of partial waves retained in the Faddeev amplitude. The convergence to 
the asymptotic behavior has also been studied. Sizeable departures are observed in 
some cases at squared momentum transfers as high as 50 (GeV/c)^. It is not clear 
whether these departures are of the order ^ or ^ log q relatively to the dominant 
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contribution and whether the bad convergence results from truncations in the cal- 
culations. From a comparison with the most complete Faddeev results, a q^ validity 
range is obtained for the calculation made in the hyper spherical formalism or in the 
Faddeev approach with the minimum number of amplitudes. 
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1 Introduction 

A great interest is currently devoted to the nucleon form factor and, in particular, to its 
behavior at high momentum transfers. Indeed, from QCD in the perturbative regime |jl|] , 
one expects the nucleon form factor, GM^q^) (or Fi(g^)), to scale like g~^, up to log terms. 
Experimentally, this behavior seems to be reached quite rapidly, around 10 (GeV/c)^ |0|. 

Many theoretical works tend to make a bridge between the low and high momentum 
domains P, ^. They deal with the non-perturbative regime of QCD, where predictions 
are more difficult, and rely on non-relativistic calculations for the lower part of the mo- 
mentum range. Among the questions that may be raised, there is the sensitivity of the 
predictions to the theoretical framework, including particular techniques, or the rapidity 
of the asymptotic behavior onset. Obviously, in view of the large momentum transfers 
that the asymptotic regime supposes, a definite statement would require a relativistic 
treatment. Some works along these lines are in progress p, ^, |^, Q. We nevertheless 
believe that a non-relativistic calculation may be of some help to provide qualitative, if 
not quantitative answers to the above questions. We will concentrate on them in the 
following. 

Quite simple descriptions of nucleons in terms of constituent quarks have relied on 
the harmonic oscillator wave function 0. In such models , the calculated form factors 
drop exponentially to zero beyond g^ = 3 (GeV/c)^ [^. Curiously, it has been sometimes 
deduced from this result that a constituent quark model could not give rise to a power law 
behavior for form factors at high q^. Better non-relativistic descriptions of the nucleon 
involve solving the Schrodinger equation, using for instance the hyperspherical harmonic 
formalism or the Faddeev equations. These approaches are often considered as exact ones, 
or almost. Calculations performed with the same quark-quark force both in Valencia 
(hyperspherical formalism) and Grenoble (Faddeev equations) |jlO[ have evidenced 
discrepancies in the binding energy of the low lying baryons of the order of a few MeV to 
be compared to a total contribution of the order of 1 GeV for the kinetic energy and the 
non-constant part of the potential [ill]. The small difference may be due to approximations 



made in either approach: the restriction to the lowest values of the grand orbital, K, in 
the hyperspherical formalism and the number of amplitudes in the Faddeev approach. In 
such conditions, it is tempting to go further and calculate the charge and magnetic form 



factors of both the proton and the neutron and see whether the discrepancy remains at 
the same level as for the binding energy. It is what we did, with the idea to check the 
sensitivity of the results to the approach and, within each of them, to the truncations that 
are currently made in the calculations. We did it also with the aim to compare the results 
to the asymptotic power law in q~^, which is expected in non-relativistic approaches with 
Coulomb or Yukawa type potentials |]T2|. While doing so, we have been led to elaborate 



simple models to understand our results. In view of their possible usefulness, some of 
them will be presented. A few remarks as for deducing the force between quarks from the 



nucleon electromagnetic form factor [|T3|, or about higher order QCD corrections will be 
made. 

The plan of the paper will be as follows. The second section is devoted to reminding 
the argument for a q~^ power law asymptotic behavior of the nucleon form factors in 
a non-relativistic approach. The origin of the difference with the g~^ QCD behavior is 
briefly explained. The third section shows the importance of the description of the wave 
function at short distances for an accurate prediction of the form factor. This is done 
on an hydrogenic type two-body system. In the fourth section, we give a few details as 
to the calculation of the nucleon wave function in the hyperspherical formalism or using 
Faddeev equations. It includes general features concerning these approaches as well as 
a few numerical results pertinent to the nucleon. Results for the form factors calculated 



with the quark-quark force of Bhaduri et al. fT^ and from different approaches are 
presented in Sect. 5. The onset of their asymptotic behavior is discussed. A detailed 
discussion about understanding some of the previous results is made in Sects. 6 and 7. 
In Sect. 8, we consider a few corrections that should be accounted for to make a realistic 
comparison with experiment. These include intrinsic quark form factor for the lower g^ 
domain of the nucleon form factors, the consideration of interaction models with some 
three-body forces, an improved description of the spin-spin force and relativistic effects 
for the higher g^ domain. As most of the work presented here deals with a non-relativistic 
picture, there is no need to introduce the extra variable Q^, which is often introduced 
in relativistic approaches to remedy the inconvenience of a possibly negative squared 4- 
momentum transfer. The following notation q"^ = (f, where q^ is also equal to the quantity 
Q^ in the Breit frame, is therefore adopted. 



2 Power law expectations for the nucleon form factor 

Predictions for the form factor of two-and three-body systems at high q have been made 



long ago by Alabiso and Schierholz |T2[ in the case of spinless constituents, assuming a 
non-relativistic as well as a relativistic treatment. For our purpose, we remind some of 
their results that may be useful for the following. 

Electron scattering on the nucleon is represented in Fig. 1. The kinematics is pertinent 
to a high momentum transfer process. In the initial nucleon at rest, quarks have a small 
momentum, essentially zero, while in the final state the three-quarks share equally the 
momentum transferred to the nucleon, q, and therefore carry the momentum |. Electron 
scattering on the nucleon involves many diagrams that differ by their time ordering and 
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Figure 1: Representation of electron scattering on a nucleon at rest (laboratory system). The 
kinematics relative to quarks indicated in the figure is that of a high momentum transfer where 
the internal momentum of quarks within the nucleon can be neglected. 



include the exchange of two gluons at least. Two of them are represented in Fig. 2. The 
first one, (a), is most often shown. The virtual photon transfers to a quark, essentially at 
rest, a momentum g that is shared with the other two quarks by the successive exchange of 
two gluons. This diagram can be considered as representing a final state interaction. The 
second diagram, (b), is completely symmetric of the first one in time. It corresponds to an 
interaction effect in the initial state and has the advantage to show that the form factor 
at high q is sensitive to the high momentum components of the nucleon wave function, a 
feature that is not so transparent on the first diagram. We insist on this point because 
interaction effects in the initial and final states have often been considered on a different 



footing in the past |15|, with the obvious idea to simplify some calculations. 
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Figure 2: Representation of two processes with gluon exchange contributions relevant to the 
nucleon form factor at high momentum transfer. The first one (a) corresponds to an interaction 
in the final state and the second (b) to an interaction in the initial state. The various -^ factors 

that contribute to the asymptotic -g- form factor in the non-relativistic approach and arise from 
the gluon or quark propagators are indicated in the figure. Many other diagrams with different 
orderings of the gluon exchanges contribute. Some of them are shown in Fig. 9. 



As expected from the previous observation, the form factor at high q is directly pro- 
portional to the high momentum component of the wave function, and not to its square. 
Let's also mention that the contribution of the two diagrams, (a) and (b), and similar ones 
not shown in Fig. 2, tend to cancel each other in an inelastic charge scattering process 
where some energy, but no momentum, is transferred to the system (this simply stems 
from the orthogonality of the states under consideration). 

Examination of Fig. 2 provides a quick estimate of the behavior of the form factor 
(or the 7*iV -^ N amplitude) at high q^. Each gluon propagator introduces a factor 
\. Furthermore, in the non-relativistic limit, each intermediate quark also introduces a 
factor \. Hence the form factor is expected to have the following behavior: 



rpn.r. 



{<f 



oc q 



:i) 



The above behavior may be invalidated if the interaction contains terms which do not 
scale like the square inverse power of the momentum at high momenta, such as a gaussian 
type force, or also if higher order effects are not small enough, in which case larger q^ 
should be considered. This is an important issue, which is aimed to be considered in 



Sects. 6, 7, 8.2 
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Figure 3: Representation in the Breit frame of a process contributing to the asymptotic form 
factor for a system of spin 1/2 particles exchanging spin 1 bosons in both a non-relativistic 
approach (nucleon case: -^ asymptotic behavior, diagram a) and a relativistic one (-^ asymptotic 
behavior, diagram b). The diagram b) only shows the extra factors responsible for a difference 
with the first case. For simplicity, the same q factor appears everywhere, but it should be 
understood that the appropriate fraction of q is to be used, depending on the diagram and on 
the position in that diagram. The representation in the lab. frame is also possible, but some 
caution is required in determining the asymptotic behavior as the time component of the four- 
momentum transfer, gO) varies like cp'. Results for the non-relativistic case also hold for spinless 
particles exchanging spinless bosons. 

The difference with the QCD power law expectation, g^"^, deserves some explanation. 
First of all, we notice that the form factor of a system of three spinless particles, interacting 
via usual scalar boson exchanges and treated relativistically, does have a behavior given 
by Eq. (1), as can be checked from Fig. 2 with the quark propagators substituted by 
scalar particle propagators and gluon exchanges substituted by scalar boson exchanges. 
However, if one considers spin 1 boson exchange, at each particle-boson-particle vertex. 



a factor q is introduced as a consequence of the vector coupling and a q~^ behavior is 
obtained. In the non-relativistic hniit, Eq. (1) holds for both types of boson exchanges. 
For quarks, spin i, the exchanges of scalar or vector bosons are also equivalent in the non- 
relativistic limit, giving rise to a q^^ power law (Fig. 3a). In the relativistic treatment, 
the spin 1 case evidences five extra factors, q, related to the definition of the quark spinors 
(see Fig. 3b), of which one has to be absorbed in the definition of the initial and final 
nucleon spinors. These factors, which arise from boosting the nucleon at rest, have a well 
defined origin and, therefore, can be accounted for, at least approximately, when a more 
realistic estimate of the form factor or a comparison with experiment is intended to be 
performed. They should not be confused with factors originating from the dynamics. The 
case of spin zero bosons coupling to spin | particles may also be considered. The above 
analysis needs to be refined but a q~^ behavior is expected too. 

Summarizing this section, a non-relativistic calculation of the nucleon form factor does 
provide a power law q~^ at high q^ (provided that the interaction between quarks is 
mediated by the usual exchange of bosons). The difference with the prediction, q~^, is 
due to both relativity and the nature of the QCD interaction, which involves the coupling 
of spin 1 bosons to spin | quarks. 



3 Asymptotic form factor of a two-body system 

It is well known that the form factor of a two-body system at high momentum transfer 
is very sensitive to its short-range description. In a few cases, a close relation can be 
established. Some results that may be found in the literature are reminded here for an 



hydrogenic type atom (see for instance [0). We emphasize points that could be relevant 



for the understanding of results concerning the three-quark system. 

Discarding the center of mass motion, the wave function of the first s-state of the 
hydrogenic atom obeys the following Schrodinger equation: 

^ ^ --e] V^(r-0 = 0, (2) 



y 2mr r 
where rrir is the reduced mass, rrir = ^ai22_ _ £ the Coulomb interaction and r = In — rol. 

' ' ' mi+m2 ' r i j- ^ i 

The normalized solution is given by: 

v47r 



while the corresponding binding energy is E = — " ™'' . Assuming now that only the 
particle 1 carries a unit charge, the form factor of this state is easily calculated: 

/.^ ^ mo 1 
df^|J\r) e''-'^^TT^ = ^— . (4) 

This expression shows that the form factor scales like q~^ at high momentum transfers. 
A direct relation of this behavior to a perturbative calculation is obtained in momentum 

6 



space. There, the form factor F{q ) gets the expression: 



-e\ / '^^ /T\ /T , ^ ""^2 



' [Ztt)'^ mi + 1712 



where (f{k) is the Fourier transform of ■?/'(r) 



^{k) = I rff e-*^^-^Xf) = v^-il^!^^^!^ = v9(-^). (6) 



(A;2 + a^m: 

— * — * 

As (^{k) is concentrated at small values of k, two domains contribute equally to the integral 
in Eq. (5), around A; = 0, and k = —q ^"^^ , which allows one to write: 

W^)-2^(g^^)/^V^(^). (7) 

mi + m,2 J (Ztv)^ 

Now, the integral in (7) is nothing but the wave function in r-space at the origin, ip{0), 
while ip{k) obeys the Schrodinger equation, which in momentum space reads: 

,P-2mrE, - /• dk' Ana ,-,, .^^ 



2mr "^' ' J {2Tif[k-k'Y 

Again using the property that (/^(A;) is concentrated at small values of k allows one to 
write: 

V{k)k-^oo = —^^(0)- (9) 

Gathering results given by Eqs. (7) and (9), one gets: 

F{q g-2^oo = 1 'ip , 10 

q^ m,2 

which is in agreement with Eq. (4) in the same limit and, at the same time, shows the 
sensitivity of the form factor to the radial wave function at the origin. This one contains 
non-perturbative effects. 

A more precise statement can be made by looking back at the expression of the form 
factor given by Eq. (4) and making an expansion of ■?/'(r) around the origin: 

^2(f) =^2(o) + 2rV'(0)^'(0) + ... (11) 

Inserting this expression in Eq. (4), one obtains from the first term a S{q} function, which 
obviously vanishes at high q. From terms beyond the first derivative, dots in Eq. (11), 
one obtains contributions that tend to zero faster than g~^ (dimensional argument). As 
a dominant contribution at high q, one is therefore left with: 

f (?),"-« = -^(^^^ii^)* «K0)«/(0). (12) 

q^ m,2 

In order to get this result, we employed the relation: 

[ / drr^ Jo{qr)]g^o = lim,^o[/ drr^jo{qr)e~''']g^o = — -. (13) 
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Equation (12) shows the sensitivity of the form factor at high q to the radial wave func- 
tion at the origin, but also to its derivative at the same point. The calculation of the form 
factor at high q therefore supposes to correctly determine the slope of the wave function 
at the origin, which, itself, is determined by the Coulomb potential, independently of the 
energy of the state under consideration. This can be checked on the first radial, / = 0, 
excitation, whose wave function and form factor are respectively given by: 

^*(f) = _L (^!!^e-^— (1 - ^^) (14) 



(l-rA2)(l-2r^^) 

In the limit q -^ oo, the above form factor exhibits for the 1/g^ term a coefficient 
different from that obtained for the ground state form factor, F{(f), given by Eq. (4), 
but this only reflects the difference in the value of the square of the wave functions at 
the origin for the states under consideration, "^^ and i°!^i:i-. For a general potential, 
the q~'^ behavior of the form factor holds provided that the interaction inserted in the 
Schrodinger equation is as singular at small distances as a Coulomb or Yukawa potential 
(oc -). It has to do with the fact that the ffist derivative of the wave function at the origin 
is determined by this piece of the interaction. It is lost when the singularity is weaker. 
It can be checked for instance that the form factor corresponding to the wave function, 
e~^^(l + br), cooked up in such a way to behave like 1 + O(r^) at small distances has no 
g~^ component at large q, the ffist non-zero component being g~®. An other interesting 
example is the gaussian wave function, ^(r) = ii^'^'^e^^ *" '^, whose form factor: 

FG^q^) = e-^'l^^'\ with 6' = 6 ^^^i±^, (16) 

1112 

has obviously no g~'^ component at large g, as expected from Eq. (12). 

While the ffist derivative of the wave function at the origin is determined by the 1/r term 
of the potential, the second one furthermore depends on the "binding energy", i.e. the 
total mass of the system minus the masses of the constituents and other constant terms 
(see Eq. (96) for the corresponding situation in the three-body case). It is therefore 
state dependent. With this respect, we should mention that hadronic systems to be 
discussed here are sensitive to a confining potential. This one does not directly influence 
the calculation of the second derivative at r = 0, but it does indirectly through its 
contribution to the total energy of the system. As will be seen in Sects. 4.3 and 6.1, 
this one is making the second derivative smaller than for a pure Coulomb problem, with 
possibly an opposite sign in some cases. This feature may have consequences for the 
rapidity of the onset of the asymptotic power law behavior of the form factor of hadronic 
systems and for the manner how it occurs. 



4 Description of the nucleon wave function in terms 
of quarks 

Determining a nucleon wave function in terms of quarks has been done in many papers, see 
refs. 0, |1|] for general presentations and refs. 0, |10|, |T|, ^, 0, H, |3|, ||, |2|, |2|, ^ 
for more specific works. We nevertheless remember a few details relative to our own 
calculations |Tl| and, especially, to some of the approximations that have been made. 



The quark-quark force employed in our calculations is that of Bhaduri et al. |ll4| : 



Vbha = o Z^ + ^ - ^ + H ^i-^i 

= E^^. = Hyin-r-), (17) 

i<j i<j 

where k = k^ = 102.67 MeVfm, a = 0.0326 MeV-^/^fm^/^ ro = 0.4545 fm, m^ = nij = 
Trig = 337 MeV, D = 913.5 MeV. With these definitions, distances are expressed in units 
of fermi and the potential in MeV. For characterizing the strength of the spin-spin force, 
we will most often refer to the following quantity: 

< = ^^ = 1.66... (18) 

Fitted on the meson spectrum, the above interaction provides a reasonably good account 



of the baryon spectrum |]TU|, although it misses the Roper resonance. It is reminded 
however that the choice of the model is not essential here and that references to other 
quark interaction models will be made in any case, see Sects. 4.3 and 8.2. 

Having neither spin-orbit nor tensor component, the force given by Eq. (17) has the 
particular feature to conserve the spin. The baryon wave function can thus be factorized 
into a spin and an orbital part. The simplification may be a drawback with some respects, 
for a discussion of the helicity conservation in QCD for instance. Here, on the contrary, it 
represents an advantage as it avoids unnecessary admixture of the effects we are looking 
at with other ones that are irrelevant for our purpose. 

Four spin-isospin quark wave functions with the spin and isospin of the nucleon, S = 
1/2, T = 1/2, are available. In notations of ref. [1^, they are: 

\s> =-^|xY + xV>, 

\A> = — =|x%^-X^r7°>, 

|M5>=i=|xV-xV>, 

|MA>=i=|xV + xV>, (19) 

where x'^iv^) ^ind x^iv^) correspond to quarks 1 and 2 coupled with an intermediate spin 
(isospin) equal respectively to and 1. All of them have a definite character under the 
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exchange of quarks 1 and 2 (symmetric : IS* > and IMS' >, antisymmetric: \A > and 
\MA > ). Under the exchange of quarks 1, 2 and 3, the spin-isospin wave functions 
IS* > and \A > are respectively symmetric and antisymmetric. The two other ones, 
IMS' > and \MA >, have a mixed character and transform into a combination of each 
other under the exchange of quarks 1 and 3, or 2 and 3, see Eqs. (127, 128). These 
wave functions have to be combined with spatial wave functions that have appropriate 
transformation properties under the exchange of quarks 1, 2 and 3 : complete symmetry 
for IS* >, complete antisymmetry for \A > and mixed symmetry for IMS' > and \MA >. 
The antisymmetry is ensured, as well known, by the color wave function, which we omit 
to write down as its effect factorizes out. The total wave function of momentum P may 
thus be written: 

"^pin, f2, fs) = iJsin, f2, rs) |S' > +^A(n, ^2, n) \A > 

+ -= (^M5(ri,r2,f3) \MS > +^MA{ri,f2,rs) \MA >) . (20) 

It obeys the Schrodinger equation: 

/ -2 -2 -2 \ 

' ^' + ^ + ^ + v{n - f2) + V{f2 - fs) + V{f3 - fi) - E] 



^ 2mq 2mq 2mq 



^p(ri,r2,f3) = 0. (21) 



For the considered potential, the center of mass motion factorizes out. Being irrelevant 
for the description of the nucleon, we omit from now on the corresponding phase factor 
as well as any reference to the total momentum, implicitly assuming that the system is 
at rest, with P = 0. As to the internal wave functions ips,A,MS,MA{'f^i-i ^^2, ''^3), they may be 
expressed in terms of Jacobi variables, fi — r2 and r^ — ^^^^^ , or any combination of them, 
see end of App. C The normalization of the nucleon wave function is chosen to be: 

d{ri - fa) rf(f3 ^~^) (V^5(^i> ^2, ^3) + V^a(^i> ^2, ^3) 

+ 2(^M5(n,r2,f3) +^MA(^l,^2,r3))) = 1. (22) 

The symmetry properties of ipMs(fi,'f^2,'r'i) and ipMA{fi,r2,r^) under the exchange of 
particles 1 and 3, or 2 and 3: 

1 v^ 

^M5(^3,^2,n) = --i^Ms{ri,f2,n) + — V^MA(n,f2,f3), 

1 v^ 

^MA(f3,f2,fi) = +-V^MA(n,^2,^3) + ^ V^M5(n, ^2, ^3), (23) 

also imply the relation: 

d{ri - f2)d{n - ^ ^ ^ ) {^Msin, f2, fs) - ^MA(n, f2, fs)) = 0. (24) 



10 



4.1 Faddeev approach 



In the Faddeev approach, the wave function ^(n, r2, f^) is written as a sum of 3 terms: 

^(ri, fs, fg) = ^12,3 + ^13,2 + ^23,1, (25) 

where \l/i2,3 is symmetric in the exchange of particles 1 and 2 while ^1^13, 2 and ^1^23,1 are 
obtained from \I/i2,3 by performing the corresponding permutations, so that '^{fi,f2,f^) 
is symmetrical in the exchange of particles 1, 2, and 3. The amplitude ^^12,3 obeys the 
equation: 



-*2 -Q 

Pi ^ P2 



-f> 

Pi 



" 2m„ 2m„ 



+ ^ + V{fi - f2) - ^)^12,3 = -Vin - f2)(^i3,2 + ^ 



32, ly 



(26) 



For calculations, \E'i2,3 is developped on the spin-isospin basis (19), or an equivalent 
one, and the spatial part, which in the present case has a total orbital angular momentum 
L = 0, is decomposed in terms of the spherical harmonics relative to the Jacobi variables 
X = n - f2 and y = -^(fs - ^^^^-j-^): 



^12,3 (n -^2,f3 



ri + r2 , 



00 



A 

E 

A=Om=-A 



y: y. yr{^)yx 



Vx{x,y). 



(27) 



The sum includes even and odd values of A. Through determining combinations with 
well defined symmetry character (see equation below), these values will be associated with 
spin-isospin wave functions that have the same parity under the exchange of particles 1 
and 2, respectively {\S >, IMS' >) and {\A >, \MA >). In practice, only a finite number is 
retained. The minimum is to retain the lowest value, A = 0, which gives two amplitudes 



l^i^i 



>, or 



in \l/i2,3 corresponding to the two spin-isospin components, \x r] > and 
equivalently \S > and IMS' >. The restriction to these two components, ipQ{x,y) and 
ipl{x, y), respectively associated with the spin-isospin wave functions Ix^r]^ > and Ix^rj^ >, 
is motivated by the expected dominance of the interaction in these channels, which should 
be somewhat corrected however for the long range of the gluon exchange. As it may be 
useful, we give here the full expression of the wave function: 



^F(n,^2,^3) 



^12,3 + ^12,3 



+ 



V2 

^^?2,3 - ^12,3 1 



1 ^ 3) + (2 ^ 3) 



+ 



V^ 



V2 2^'-'^-2^'-'^ 

"^32,1 '~ "^32,1 /, 



S> 

\MS> 



V2 



\MA>, 



:28i 



where, in the two component case: 



;0,1 



,0,1 
V'13,2 



0,1, 

V'o ( 



0,1/ 
^0 ( 



\-x H y\ 

'2 2 ^' 



1 



-y H x\ 

2^ 2 ' 



-X 



Vs 



y\ A -nV- 



Vs. 



-x\ 



(29) 
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While values of A higher than are neglected in the above minimal calculation of the 
Faddeev amplitude \E'i2,3, their contribution to the total wave function are not totally 
neglected. They appear through the exchange terms \E'32,i and \l'i3,2, whose expansion 
in terms of spherical harmonics J2Yy'*{£) ^y iv) ^-nd spin-isospin wave functions, in 
principle infinite, has been limited to 10 terms in calculating the observables (form factors) 
presented in this paper. A calculation consisting in retaining 8 amplitudes in the expansion 
of the Faddeev amplitude \l'i2,3 (not to be confused with the wave function) has also been 
performed. The comparison with the 2 amplitude calculation will provide an interesting 
test. Indeed an important issue is to know whether some truncation in the Faddeev 
amplitude affects the power law behavior of the form factor at high momentum transfers. 



4.2 Hyperspherical formalisin approach 

In the hyperspherical formalism, a different choice of the variables is made (see App. A). 
The modulus of the Jacobi variables are expressed in terms of an hyperradius, p, and an 
extra angle, : 

X 

-j= = Pi2 = P sin0, 

— = = p3 = p COS0. (30) 

In terms of these variables, the total wave function reads: 

M/HH = l^(p,fi) >= E ^kAp) K'',S(^) l^y^ >]j, . (31) 

K,L,sym 

where \sym > represents the spin-isospin wave functions given by (19) and Q the set of 
the various angles relative to the unit vectors, pu and ps together with the hyperspherical 
angle 0. ^fx!s«mi(^) stands for the hyperspherical harmonic (HH) of definite symmetry, 
orbital angular momentum and quantum numbers specified by K. The hyperradial wave 
functions, 4'k,l{p), satisfy a (infinite) set of coupled equations written in detail in App. 
A. 

In practice, as for solving the Faddeev equation, a truncation is required. Two ampli- 
tudes have been retained. They involve the totally symmetric spin-isospin wave function 
IS* > and a combination of the mixed ones, IMS' > and \MA >. They respectively have 
the lowest allowed K value, and 2. By identifying 

i^o,L=oip)=Mp), (32) 

^2,L=o(p) = ^3(p), (33) 

an explicit expression of the wave function is: 
*HH(ri,f2,f3) ~ —^ [Mp)\S > 



+72 ^3(p) ((cos^0 - sm'^(j))\MS > -2 sin0 cos0 p^.pslMA >)] . (34) 
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The front factor ensures the normahzation over the various angles 



^3 



1, (35) 

and the radial wave functions are normalized as: 

dpp'{ij!{p)+r3{p)) = l. (36) 



As for the Faddeev approach, one may wonder whether the truncation of the wave 
function has some consequence for the prediction of the form factor at high momentum 
transfers. From now on, it is remarked that the description of the system with the 
totally symmetric radial wave function ipi{p) is more economical but poorer than the 
corresponding one in the Faddeev approach, first bracket on the r.h.s. of Eq. (28). The 
latter contains extra terms that, in the hyperspherical formalism, suppose to introduce 

components with K = 4, 6 This can be seen by realizing that the completely symmetric 

hyperspherical harmonic for K = only contains Yx=o{x) Yx=o{y), whereas the Faddeev 
amplitudes \E'32,i and \E'i3,2 in the simplest case (2 amplitude calculation) contain terms 
with A 7^ 0. Let's also notice that the radially symmetric component, ipiip), is the 
dominant one, while the mixed symmetry component, ipsip), which is smaller, is directly 
determined by the spin-spin interaction. 



4.3 Static properties and wave functions 

We give in Table 1 a few results relative to the nucleon wave function calculated with 
the Bhaduri et al.'s force in the different approaches : Faddeev with 2 and 8 amplitudes 
and hyperspherical formalism. They concern the mass of the nucleon and its square 
matter radius, the proton and neutron charge squared radii and the mixed symmetry 
probability. Results where the spin-spin force in the hyperspherical formalism is neglected 
are also given. The comparison of the results incorporating the spin-spin interaction 
does not evidence much difference for the mass or matter radius. The sensitivity to the 
approximation is more important for the neutron charge squared radius and the mixed 
symmetry probability, which both involve the spin-spin interaction at the first and second 
order respectively. 

As we are interested in the form factor at high q, which depends on the description 
of the system at short distances, we also show some wave functions. This is done only 
for the hyperspherical formalism where the graphical presentation of the results is much 
simpler than for the Faddeev approach since it reduces to the wave functions ipiip) and 
ipsip). They are given in Fig. 4. Examination of V'i(p) at short distances evidences at 
first sight an exponential behavior. This one is approximately given by: 

^,{p)=a,e-^^P, (37) 

with ai = 19.95 fm"^ and Pi = 1.715 fm"\ 

The slope, f3i, which is entirely determined by the Coulomb like part of the potential 
(including the spin-spin part which also behaves as - at small r), is close to what is 
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Table 1: Static properties pertinent to the description of the nucleon, as calculated from different 
approaches with the same force: Faddeev equations with 2 and 8 amplitudes (first and second 
columns, F) and hyperspherical harmonic formalism (HH) with a restriction to the K values, 
and 2 (third column). Results calculated in the hyperspherical harmonic formalism with 
K = Q are given separately in the fourth column for a different force. This one only includes 
the Coulomb and confining parts of the Bhaduri et al.'s potential (no spin-spin component). 
The listed quantities are successively the mass, the square matter radius, the mixed symmetry 
probability, the proton and neutron square "charge" radius, and the density of the nucleon at 
the origin. 

(F) (F) (^) pfflT 



M;v(MeV) 


1031 


1020 


1039 


1201 


<rl> (fm^) 


0.219 


0.219 


0.218 


0.255 


Ps' 


1.0% 


2.1% 


1.5% 





< r-l > (fm^) 


0.234 


0.243 


0.238 


0.255 


< rl > (W) 


-0.015 


-0.024 


-0.020 





^2(0)(fm-6) 


386 


386 


398 


136 



theoretically expected (see App. A): 

8^2 



5n 



K + K'^)mg = 1.702im-\ (3^ 



Some departure from the exponential form expected in a pure Coulombian problem 
is however observed, which is seen in particular on the wave function, ipi{p), calculated 
when the spin-spin part of the interaction is turned off (see Fig. 4). For this case, a 
sizeable change occurs in the slope of the wave function at short distances when going 
from p = to p = 0.25 fm. It can lead to specific features in the form factor at high q, 
which, as reminded in Sect. 3 for the two-body case, is sensitive to that part of the wave 
function. Numerically, this one up to p^ terms is given by: 

Mp) = (^'iil-P'iP + l'ip'), (39) 

with a[ = n.65im~^,P[ = 0.645 fm"^ and -f[ = -2.0fm~l 

The value of P[ is close to the theoretical expectation, ^^nniq = 0.640 fm~^. The value 
of a'l is smaller than the value of ai in (37). This is in direct relation with the absence of 
the short range attraction that the spin-spin interaction provides on the average. As to 
the '-/[ coefficient, it evidences a striking departure from its value in the pure Coulombian 
problem, 7 = ^ = 0.20 fm~ (according to the notation employed in App. A), since it 
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,)(rm) 



p(rm) 



Figure 4: Radial wave functions V'i(p) ^^'^ V'3(p) obtained in the hyperspherical formalism with 
the Bhaduri et al.'s force, Eq. (17) (continuous line). The wave function ipi{p) in the case where 
the spin-spin force is turned off is also shown (dashed line) . 

has an opposite sign and is one order of magnitude larger. It depends on the constant 
terms appearing in Eq. (88) and should be compared to the theoretical expectation given 
by Eq. (96): 



-^'- 



3m 



,+2^: 



5 / 



3^2 



-Kiria 



-2.1 fm" 



(40) 



Through E, it mainly involves the non-perturbative part of the potential (17), arising 
from the confinement. In comparison, the role played by the perturbative part due to 
one-gluon exchange, represented by the last term in the l.h.s. of (40), is negligible. 

The behavior of the component, ip3{p), at short distances also deserves some comments. 
An approximate expression is: 



'ipsip) = Oispe 



-PsP 



(41) 



with ^3 = 8.4 fm , P^ = 3.76 fm . The value of a^ may be compared to the theoretical 
expectation given by Eq. (101): 



32 

35^ 



K^ruqai 



8.23 fm~ 



(42) 



As to the (3^ parameter, it has no well determined interpretation. While the exponential 
factor in Eq. (41) turns out to be a (very!) good approximation to the numerical solution 
at small p, a more complete theoretical analysis indicates that the second term of the 
expansion of the solution should contain a well defined term, p^ log p, besides a p^ term. 
While the former is determined by the knowledge of 0:3, ai and /9i, the latter is related 
to the solution of the homogenous Eq. (100) and cannot be similarly determined. The 
insertion of the p^ term in the l.h.s. of this equation indeed leads to an indertermination 
for its coefficient. 
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The comparison of the phenomeno logical parameters f3i, f3[, '-/[ and 0:3 with the theo- 
retical expectations evidences some slight discrepancy. We believe that this one has its 
origin in the numerical methods used to solve Eqs. (99, 100). Notwithstanding it, we 
think that the results show the adequacy of the Numerov algorithm that we employed 
to solve the equations of the HH approach and determine the short range part of the 
hyperradial wave function. 



Expression and calculation of the nucleon form fac- 
tors 



Having determined the nucleon wave functions, we now consider the form factors that 
can be calculated from them. The charge and magnetic form factors of the proton and 
the neutron, which represent four quantities, are the object of a current interest. They 
can be obtained from the charge and magnetization densities. In the approximation used 
for describing the wave function, where the total orbital momentum of the nucleon is 
L = 0, in agreement with the quark-quark force we used, the starting point is given by 
the following formulas: 

< m^^GW) + ^-^GUcf)\N > {2T,f5{pf -P'-q) = 

< N{P^)\j:ilil + r/)e^«-)|iV(P^) >, (43) 

< iV|a(i±^^Ml!l + l^l!^M^)|iV > (27r)3<5(p/ - P' - q) = 

' ^ 2 2mjv 2 2mjv ^1 ^ ^ ^ ^^ 

< N{P^)\i: ^|(^ + r/)e'^-|iV(PO > . (44) 

Using the symmetry of the nucleon wave function, the sum over the three-quarks ap- 
pearing in (43, 44) can be restricted to one of them, the matrix element being multiplied 
by 3. The third quark is here chosen as the expression of our wave functions for mixed 
symmetry states have a definite symmetry property under the exchange of particles 1 and 
2. After performing the algebra relative to spin and isospin operators, one is left with the 
following expressions: 

G-P (g®) = (< S\0\S >+- < MS\0\MS >) + V2< S\0\MS > 

+ (< A\0\A >+-< MA\0\MA >) - V2 < A\0\MA >, (45) 

Gl{q^) = -V2<S\0\MS>+V2<A\0\MA>, (46) 

Gliiq^) ^ = {< S\0\S + l < MS\0\MS >) + V2 < S\0\MS > 

TTljq 2 

--{< A\0\A >+- < MA\0\MA >) + — < A\0\MA >, (47) 

m 2 1 \/2 

Guif) -^ = — (< SIOIS > +- < MSIOIMS >) - — < SIOIMS > 
mjv 3 2 3 

16 



+^(< A\0\A >+\< MA\0\MA >) + ^ < A\0\MA > . (45 

The matrix elements are defined as: 



<X\0\Y>= jdC^^^)d{r,- 



2 ' 
^x(ri, f2, fg) e*^"-(3^^3-^V^) ^y(fi, f2, fg), (49) 

where X, F stand for the different symmetry character of the components of the wave 
function, S, A, MS and MA. The expressions (45-48) have been written in such a way 
to emphasize their dependence on four independent quantities. As one of them involves 
the completely antisymmetric component ipA{'^i,'^2,'f^3) which is expected to be small 
(when it is not simply put to zero as a result of approximating the wave function as 
in the hyperspherical formalism or the two amplitude Faddeev calculations), one should 
approximately obtain the following relation: 

\i GW) - Gli>f) ) - G^^icf) ^ - G^m ^ = 0. (50) 

3 mjv rriN 

This is trivially satisfied for: 

GW) = fief), Glif) = 0, 
G%if)^ = /(a, Gl,{q^)^ = -Ifif), (51) 

a choice which is sometimes made in the litterature, f{(f) being given most often a dipole 
expression. The validity of Eq. (50) has been tested in ref. [^ for both the Faddeev 
wave functions used in the present work and the measurements. At g^ = 0, expressions 
different from Eqs. (45-48) may be found in the literature 0. They correspond to another 
type of component \A > with a non-zero orbital angular momentum (1 instead of 0). In 
the harmonic oscillator model, the two kinds of states represent at least a 2huj {L = 1) 
and a 6huj [L = 0) excitations respectively. 

The results for the form factors, G^^{(f), G]^{(f), G^ji^cf) and G\j{(f) are given in Fig. 
5, for transferred momenta up to g^ = 5 (GeV/c)^. In each case, the results are shown 
for the Faddeev calculations with two and eight amplitudes and for the hyperspherical 
formalism (with the two waves i^ = and K = 2). For G''^{(f), we also give results 
for a gaussian wave function corresponding to a matter radius close to that obtained 
in the other approaches. It is seen that the Faddeev calculations with two and eight 
amplitudes significantly differ from each other, especially for the charge neutron form 
factor. The eight amplitude results compare reasonably well with the results obtained in 
the hyperspherical formalism up to g^ = 3 (GeV/c)^. Some discrepancy appears beyond 
this value, as seen in this figure. One also sees that the result using a gaussian wave 
function quickly differs from the above ones (a factor of about 2 at g^ = 5 (GeV/c)^ ) and 
tends to more rapidly. 
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Figure 5: Electric and magnetic form factors of the proton and the neutron up to g^ = 
5(GeV/c)^. Except for G^{(i^), where it is not necessary, they have been multiphed by a 
factor q'^ to emphasize the differences in the higher q domain. Results are presented for wave 
functions calculated with the Faddeev equations (2 and 8 amplitudes) and the hyperspherical 
formalism (K=0 and 2). A comparison with results provided by a gaussian wave function can 
be made by looking at the upper-left figure (dotted line). 
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Figure 6: Same as in Fig. 5, but up to g^ = 50(GeV/c)^. Results liave been multiplied by q^ 
in all cases to emphasize the asymptotic behavior. 
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In Fig. 6, we present results that concern the asymptotic behavior of the form factors. 
They have been divided by q~^, which is the expected behavior at high q (Sect. 2). These 
results, which should tend to a constant number, are shown up to g^ = 50(GeV/c)^. 
Examination of the figures indeed indicates that, roughly, some asymptotic value seems 
to be reached. This demonstrates that constituent quark models, contrary to what is 
sometimes claimed, do lead to a power law behavior of form factors. This is achieved pro- 
vided that they incorporate at least a minimal description of the short range correlations 
produced by gluon exchanges. Looking in more details, one may however notice that the 
results obtained with the hyperspherical formalism are somewhat higher and that their 
convergence to some asymptotic value is less clear than for the Faddeev approach. In this 
case, a slight variation with (f is still visible, suggesting that the asymptotic behavior may 
not be reached yet. These different features have led us to make a more refined analysis. 

As a side remark, let us mention that the smallness of the form factors at the highest 
values of g^ under consideration (of the order of 10~^ — 10~^) may cast doubt about their 
numerical accuracy. Tests for which an analytical result was available (see next section) 
has revealed that this accuracy was surprisingly good and could not be responsible for 
the above departures. In fact from a dimensional argument, the form factor given by Eq. 
(49) is expected to at least contain a factor, -^ (see App. B), which already explains for a 
large part the smallness of the form factor. Some of the observed departures may however 
be due to the lack of accuracy in solving equations to obtain wave functions. 



6 Discussion of the results 



Due to the difficulty to perform extensive calculations with the Faddeev approach, the 
discussion made below will essentially concerns the results obtained with the hyperspher- 
ical harmonic formalism. We will nevertheless attempt to establish some relationship 
between the two approaches. On the other hand, instead of working with the proton and 
neutron, charge and magnetic form factors, we will deal with matrix elements involving 
components of the wave function with a given symmetry, namely: 

<5|0|5>, <S\0\MS>, <MS\0\MS>, <MA\0\MA>, (52) 

whose expressions are given by Eq. (49). The reason to make this alternative presentation 
of the results is the possibility that they may not all have the same asymptotic behavior. 
This is expected from other three-body calculations where it was shown that, in momen- 
tum space, the component ipMS was going to zero more rapidly than ips HH- This result 



was however obtained with separable forces and, actually, for the Faddeev amplitude. 



6.1 Calculations with the hyperspherical harmonics 

The advantage of the hyperspherical harmonic formalism is to provide a simple, but 
approximate expression of the wave function at short distances, which furthermore can 
be used for analytic calculations of the form factors. Thus, from the expressions ipiip) 
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and ipsip) given by Eqs. (37) and (41), one gets for the various matrix elements listed 
above: 



<^|0|^> = 405^6 



al(3i 



(g2 + 6^2)7/2 



'•'^GeV/c)^ 



<S\0\MS> = -2835^6- ^^ 



13.5 



(g2 + 6( 



hf33\2 



)9/2 



On 



< MS\0\MS > = 7290^ 64 



Oq 



< MA\0\MA >= -2430^ 64 - 



(g2 + 6/51)11/2 

168V6(3sq''> + .. 
(g2 + 6/51)11/2 



q 

75.6 
25.2 



g« 



f(GeV/c)^ 
(GeV/c)^ 
(GeV/c)^ 



(53) 
(54) 

(55) 
(56) 



For the two last matrix elements, only the terms dominant at high q and the first cor- 
rections to it have been retained here. A complete expression is given in App. B. An 
expression identical to Eq. (53) was obtained in p3|, Q, but on a pure phenomenological 
basis, without relation to some microscopic dynamics. 





q <MA|0,MA> 



q" (GevVc') 

a) 



q' (GevVc') 

b) 



Figure 7: Form factors calculated in the hyper spherical formalism up to q'^ = 50 (GeV/c)^. The 
form factors, < S'|0|S' > and < S\0\MS >, have been multiplied by g^ (Fig. a) and the form 
factors, < MS\0\MS > and < MA\0\MA >, by a factor q^ (Fig. b). In each figure, we present 
the form factors obtained from the numerical wave functions, ipi{p) and ipsip) (continuous line), 
as well as the analytic ones, Eqs. (53-56) (dashed line). The < 510 IS" > form factor in absence 
of spin-spin interaction is also shown in Fig. a (dash-dotted line for the numerical calculation 
and dashed line for the analytical one). 

Being a good representation of the short range behavior of the calculated wave functions, 
the parametrizations of the components, V'i(p) ^"^^ '4'z{p)i given by Eqs. (37) and (41) 
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should provide accurate expectations for the form factor at high q in the hyperspherical 
harmonic approach considered here. By comparing the form factors given by Eqs. (53-56) 
to those calculated with the actual wave functions, it is possible to determine the range 
of (f where the description of the wave functions by the parametrizated ones, Eqs. (37, 
41), becomes relevant, as well as the role of non-perturbative effects they do not account 
for. The results are presented in Figs. 7a, b. They have been multiplied by a factor g^ for 

< S\0\S > and < S\0\MS > and q^ for < MS\0\MS > and < MA\0\MA >. These 
factors are chosen in accordance with the asymptotic behavior expected from Eqs. (53- 
56). The results so obtained can be directly compared to the asymptotic numerical values 
given in Eqs. (53-56). Figure 7a also contains the result for the form factor, < 5*1 OIS" >, 
when the spin-spin part of the Bhaduri et al.'s force is removed. 

Examination of the figures shows a good agreement between the form factors given by 
Eqs. (53-56) and those obtained with the calculated wave functions, ipi{p) and i'sip), in 
the range, g^ > 30 (GeV/c)^. With the help of the analytic expression, it is possible to 
have a discussion on the onset of the asymptotic behavior. 

These results evidence two striking features. First, the form factors < 5*1 015* > and 

< S\0\MS > at very high q behave hke g~^, instead of q~^ as it is expected in a complete 
calculation [1^. The g~^ behavior can be understood using counting rules, similarly to the 



derivation of the q~ in the general case: the propagator of the set of three-quarks in the 
hyperspherical harmonic formalism introduces a factor g~^, while the interaction, given 
by the Fourier transform of the - hyperradial potential provides a factor q~^. Second, the 
form factors < MS\0\MS > and < MA\0\MA > behave like q~^, but the asymptotic 
behavior is far to be reached in the range around 50 (GeV/c)^ (a factor 2-3 is missing). 

The explanation for the wrong behavior of the form factors < S'lOIS' > and < 
S\0\MS > is to be looked for in the approximation consisting in retaining the lowest 
K values in the expansion of the wave function. Let us examine for example the spatial 



wave function for K = 0, V'i(p). Depending on the variable p = J Pu + pi, at small pu 
(Pi2 << Ps ) it behaves like : 

Mp)p^.^o = MPs) + IMps)— + •••• (57) 

^ Ps 

The appearance of PI2 as the first non zero term in the expansion cannot account properly 
for the short range correlations of particles 1 and 2, the third particle being a spectator. 
As one obtains in the two-body case (see Sect. 3), a term linear in p^ is expected. To 
get the right behavior, one should sum up an infinite set of contributions with K ranging 
from to 00. Furthermore, a rapid convergence is not ensured at all. 

In the case without spin-spin force [ipsip) = 0), we looked at some contributions of 
the next terms in the hyperspherical harmonic expansion corresponding to i^ = 4 and 
6, Eq. (102)). The first one contributes to < S'lOIS* > three times as much as the wave 
i^ = in the asymptotic domain, with the same sign. The contribution of the K = Q 
wave, which is antisymmetric in the pi2 and ps variables and, thus, introduces a richer 
structure in the wave function, could correct the result in the right direction. This is 
in agreement with the fact that reproducing the correct asymptotic behavior requires a 
somewhat fine description of the wave function at short distances, especially with respect 
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to the hyperspherical angle. 

Another related aspect is worthwhile to be mentioned. Up to a factor, the coefficients 
of the q~'^ term in Eq. (53) and Eq. (54), ai.aiPi and ai.a^ respectively, have a form 
similar to Eq. (12), namely the product of the wave function at the origin, ai, which 
by itself contains non-perturbative effects, and its derivative at the same point, aif3i or 
0:3. Factorizing out the wave function at the origin, ai, the last quantities, j3i and ^, 
for which an analytic expression can be obtained from Eqs. (38) and (42), are seen to 
be linearly dependent on the strength of the dominant part of the interaction at short 
distances, determined by k, and k,'^. This is not in accordance with the argument reminded 
in Sect. 2, which implies two gluon exchange and, therefore, a dependence on the square 
of these k, factors. 

On the other hand, from the smallness of the mixed symmetry state (1.5%), it may 
be inferred that its contribution to the form factor is negligible. This is true at small 
momentum transfers, but not at high ones, as shown in Fig. 7a or on the analytic 
expressions given by Eqs. (53, 54). Taking into account the expressions of a^/ai, Eq. 
(42) and /3i, Eq. (38), it is found that the matrix elements < S\0\MS > and < S'|0|S' > 
at high q are quite comparable. Indeed the ratio tends to a constant of the order of unity: 

^ <s\o\s> ^^-- - -^^2 ;7T< " '^-^ ^^^^ 

The last aspect we want to discuss about the form factors, < 5*1 OIS* > and 

< S\0\MS >, is the role of non perturbative effects. In a pure Coulombian problem 
and in the approximation where only the K = wave is retained, the form factor ap- 
proaches its asymptotic value from below. Examination of Fig. 7a shows the opposite 
feature, as a consequence of having incorporated into the p^ term of the description of the 
wave function at short distances a contribution depending on the energy of the system, 
see Eq. (40), which indirectly involves the confining interaction. The above feature, which 
is strongly enhanced when the spin-spin force is neglected (see Fig. 7a), will be also seen 
in other calculations presented below. The observed structure at g^ = 4 (GeV/c)^ in the 
case without spin-spin force has to do with the change in the second derivative of the 
wave function around p = 0.3 fm that can be seen in Fig. 4. An immediate consequence 
is to make the onset of the asymptotic behavior to be reached quicker. 

It is also interesting to look at the two other form factors, < MS\0\MS > and < 
MA\0\MA >, although it is more difficult to rely on theoretical expectations, due to 
the lack of interpretation for the parameter (3^ entering the short-range description of 
the component ipsip), Eq. (41), on which they depend. We already noticed that these 
form factors behave like q~^ at high q (see analytical calculations, Eqs. (55, 56)). On 
this basis, one may think that their role should be negligible at high q, in comparison to 

< S\0\S > or < S\0\MS >. This is hardly verified (< MS\0\MS >= 0.55 < S\0\S > 
at g2 = 50(GeV/c)2). 

Moreover, the values of these form factors multiplied by q^, which should tend to a 
constant, are still increasing in the range g^ ~ 50 (GeV/c)^ and far to have reached the 
asymptotic value (larger by a factor 2). 
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In fact, assuming that the asymptotic behavior was reached in this range led us to the 
prehminary conclusion that this form factor was behaving like g"^, like < S'|0|S' > and 

< S\0\MS >. It is only from the examination of the analytical calculation, Eqs. (53-56), 
that we could determine that the asymptotic behavior of the form factors, < MS\0\MS > 
and < MA\0\MA >, was q~^ and understand, as well, the other features they evidence. 
In particular, as is seen from Eqs. (55, 56), the next to leading order contribution in 
g~^ contains an extra factor q~^, instead of q~^ for the form factors, < S'|0|S' > or 

< S\0\MS >. The contribution is about 60% of the dominant term at g = 7 (GeV/c) 
and, since it is destructive, it provides an overall reduction of the form factor by a factor 
2-3. This feature strongly delays the onset of the asymptotic behavior for the form factors, 

< MS\0\MS > and < MA\0\MA >. Equations (55, 56) suggest that the asymptotic 
behavior is obtained within a few % at values of g^ as large as 10000 (GeV/c)^ (!). 

So, we can conclude that calculating reliably the form factors at values of g^ as high as 
those considered here supposes not only an accurate determination of the wave function at 
short distances, but also on the whole range. Calculations we did up to g^ = 200 (GeV/c)^ 
have revealed moderate but regular oscillations around the asymptotic behavior of the 
form factor < 5*1 OIS* >, explaining a slight plateau that can hardly be seen in Fig. 7a 
around g^ = 40 (GeV/c)^, but shows up in the numerical values. From the checks we 
did, it turns out that these oscillations are in relation with the value of the hyperradius 
where the inner and outer Numerov solutions of the Schrodinger equations were matched. 
At this point, the third derivative of the numerical solution evidences a change in sign 
(giving rise to the oscillations) that does not seem to have any physical meaning. 



6.2 Calculations in the Faddeev approach 

One may wonder whether observations made about the calculations performed in the 
hyperspherical harmonic approach apply to those in the Faddeev approach, for which we 
have a priori no precise benchmark as far as the asymptotic behavior of form factors 
is concerned. From examining the corresponding results presented in Fig. 8, we may 
notice that the convergence of the form factors < 5*1 OIS* > and < 5*10 IMS' > to some 
asymptotic behavior assumed to be q~^ is the same as for the hyperspherical harmonic 
calculations with respect to g"^: same rate and from above for the first one. As to 
the form factors < MS\0\MS > and < MA\0\MA > (+2 < A\0\A >), the absolute 
value of their product by g^, especially for the second one, tends to steadily increase 
around g^ = 50 (GeV/c)^, quite similarly to the results obtained with the hyperspherical 
harmonic approach. Due to the way calculations were performed, the contribution of the 
fully antisymmetric spin-isospin state to the last matrix element, < A|0|A >, could not 
be easily separated in the case of the 8 amplitude Faddeev approach where it is not zero. 
Its contribution is in any case quite small and should not affect the discussion relative to 
the matrix element, < MA\0\MA >, to which it is admixed (what we reminded between 
parentheses) . 

Two more observations are in order. There is no evidence that, for g ^ oo, the form 
factors, < MS\0\MS > or < MA\0\MA >, go to zero more rapidly than the form 
factors, < S\0\S > or < 5*10 IMS' >, as is expected from some calculations (see Sect. 7.1 
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Figure 8: Form factors calculated in the Faddeev approach up to q^ = 50 (GeV/c)^. The dashed 
and continuous lines respectively represent the results for the 2 and 8 amplitude calculations. 
The form factors, < S\0\S >, < S\0\MS >, < MS\0\MS > and < MA\0\MA >, have been 
multiplied by a factor q^. The values, which these quantities are expected to converge to, are 
given in Eqs. (74). 

below). On the other hand, on the basis of the same calculations, the form factors in the 
Faddeev approach are expected to go to zero more rapidly and to be smaller in absolute 
value at high q than in the hyperspherical harmonic approach with the truncation space 
chosen. The comparison of the results indicates that this conclusion is supported in some 
part for the form factors, < ^|0|^ > and < S\0\MS >, (a factor 0.4 at q'^ = 50 (GeV/c)^ 
for < S\0\S > for instance), but not for < MA\0\MA > (+2 < A\0\A >) (a factor 2.5). 
At this point, the slow convergence of this matrix element to an asymptotic value prevents 
one to make definite statements. 

Another question of interest concerns the convergence of the form factors to their asymp- 
totic behavior depending on the Faddeev calculation, two or eight amplitudes in the 
present case. Apart for a change in scale for the form factors involving the mixed sym- 
metry state, the examination of our results does not evidence any significant change in 
the asymptotic behavior. The increase with q of the quantity, q^.{< MA\0\MA > (+2 < 
A|0|A >) ), is slightly less with the eight amplitude than with the two amplitude Faddeev 
calculation around g^ = 50 (GeV/c)^. In any case, and contrary to the hyperspherical har- 
monic approach, there is no evidence at this point that the power law asymptotic behavior 
of form factors is affected by the truncation of the Faddeev type calculation which, in the 
simplest case (two amplitudes), already involves a more structured wave function. Only 
the overall size may depend on the truncation. In fact, an examination of the mixed 
symmetry components of the wave function at small distances for the 2 and 8 amplitude 
cases, around p = 0.05 fm, shows a discrepancy of a few percents only, without any re- 
lation to what is observed in the calculations performed for momentum transfers up to 
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q^ = 50 (GeV/c)^. This is another indication that the asymptotic behavior is far to be 
reached for the matrix elements < MS\0\MS > and < MA\0\MA >. 



Wave function approach to the asymptotic form 
factors 



Partly to reduce the gap between the two approaches we discussed, partly to answer some 
of the questions raised by the results obtained in the Faddeev approach, we here present 
further developments. The aim is to get an estimate of the asymptotic form factor as 
accurate as possible, starting from the observation that the wave function at the origin is 
rather well determined and independent on the approach (see Table 1). Such a program 
can be in principle performed because, similarly to the two-body case discussed in Sect. 
3, the form factor at high q is expected to only depend on the gluon exchange force, once 
the wave function at the origin is known. The expression for the form factors reads in 
momentum space: 

< ^|0|r >= / ^^ ^x(^i2, /?3) M^i2, '^a + \/| g). (59) 

Being essentially interested in the asymptotic behavior, some approximation can be done. 



We follow the qualitative lines developped by Alabiso and Schierholz ||12|, whose ap- 
proach essentially extends to the three-body system that one reminded for the two-body 
case in Sect. 3. Noting that the wave function is peaked at small values of the argu- 



ment, two regions of the integration range over the variable K3 are important: K3 = and 



K3 = —\\q. One thus gets : 



< X\0\Y >,_= I / i|ji|2j *v(S,„S3)) ^V(0, ^l ,-) 

dR^o dK,r, ,^ ^ A , /2 ^ , , 

V'y('^i2,/?3) ^x{0,-\-q). (60) 




(27r)3(27r)3 ^^ ^ ^"' "'J ^"^ ' VS 

Remembering that the integral is nothing but the value of the wave function in con- 
figuration space at the origin, and that only ips{p = 0) 7^ 0, we come to the conclusion 
that two among the four matrix elements < 5*10 15* >, < S\0\MS >, < MS\0\MS > 
and < MA\0\MA > have a non zero contribution at the dominant order, which is given 

by ipY{i^i2 = 0, /?3 = yi g) or ipxii^u = 0, K3 = -•/§?)• They are: 

< S\0\S >,^oo= 24\p = 0) 4')(/5i2 = 0, /?3 = ^Jlq), (61) 

< S\0\MS >,^oo= 4\P = 0) ^Sk'^is = 0, /?3 = y^g). (62) 

The superscripts, r and k, remind that the corresponding wave functions, which have 
been defined previously, refer to different spaces. 
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7.1 Improving upon the hyperspherical harmonics results 

One can check that when Eqs. (61) and (62) are apphed to the numerical solutions of the 
Bhaduri potential at short distances obtained in the hypersperical harmonic formalism, 
Eqs. (37) and (41), the asymptotic behavior of the form factors, Eqs. (53) and (54), is 
reproduced. Concerning < MS\0\MS > and < MA\0\MA >, the next to dominant 
order approach confirms the asymptotic behavior obtained in Eqs. (55) and (56). 

In order to improve the calculation of the form factors in the hyperspherical approach, 
we can try to get a better determination of the wave function for high values of the 
momenta, what can be done iteratively from a zeroth order wave function. This one is 
chosen as the Fourier transform of a wave function of the form ■?/'i(p)/7r^/^ = ae~^'' /n'^^'^ 
(it may correspond for instance to a Coulomb like potential or to the solution of the 
Bhaduri et al. potential at short distances Eq. (37), in which case, a = ai, (3 = (3i). It 
reads: 



^W(;5,2, /?3) = j dpi2 dps e-'^^^'-^^'^''''-^'^ -^^ e-'^V^?^+^ 

iVV2 



('^f2 + «i + ^')'/'' 



(63) 



where 



/tl2 



^^■-P= s, . . I (ft _ ^±A), N^n = n,.^n-^-fs. 



•n 

The quantities, pi,p2 and ps, represent the momenta of the quarks 1, 2 and 3. Re- 
presenting only the high momentum behavior of the wave function, '?/'*-°-'(/?i2, ^^s) and 
subsequent wave functions don't fulfill the usual normalization condition, which would 
imply a relation between a and /^, see Eq. (98). 

For an interaction dominated by a Coulomb force (or a Yukawa one) and for the wave 
function given by Eq. (63), it is possible to calculate analytically the wave function by 
an iteration procedure up to first order. Thus from the Schrodinger equation, 

we can write (as it is done in the Faddeev formalism): 

^ = V'12 + V'13+V'23, (65) 

SO that 

^ = ^i^V^^^^'^- (66) 

The contribution due to the force between particles 1 and 2, which can be identified to 
the Faddeev amplitude V'12,3, Eq- (27), but in momentum space, reads: 



X 



'12, ^^3) 


2v^. 
- 15. ^' 

3iVi/2 

■12 + ^^3 + 


+ <) 

^2)7/2 


rrig 
log( 
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+ /tl2 
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3 iVV2 N^ \ 

+2 {Kh + 4+ p'n>^i + p') ^ {^2 + 4+ P'n4 + p')v ■ ^^ ' 

The total wave function is obtained by adding to (67) the contributions where the role 
of particles 1 and 3 or 2 and 3 are exchanged, corresponding to the interaction of the pairs 
23 or 13. As a consistency check, it can be verified that averaging the total wave function 
over the hyperspherical angle with the appropriate weight factor in the differential 
element of volume, sin^0 cos^0, allows one to recover the zeroth order wave function (63). 
This is achieved provided that the following condition is fulfilled for a Coulomb potential 
with strength, k, + k'^: 

8^2 
bn 



P = ^iK + K'^)mg. (6^ 



Expression (67) can also be used to calculate the mixed symmetry component. In such 
a case, the separate amplitudes with the pair of particles 12 having spin and 1 should 
be associated with the appropriate strength of the force, k + 3k,'^ and k — k'^. The Yukawa 
nature of the spin-spin force may also be accounted for. In our opinion, the wave function 
so obtained may be usefully compared to the two amplitude Faddeev calculation presented 
in Sect. 4.1. It certainly has not all the physics contained in this one, especially that 
related to the confinement but it probably has a large part of the physics relative to the 
description of the form factor at high momentum transfers. In this respect, it is noticed 
that the amplitude (67) behaves like k^ when k^ goes to oo, while the zeroth order 
amplitude (63) behaves like k^ . As it will be shown below, this property determines the 
behavior of the form factor at high q for the dominant contributions. 

The expression of the total wave function (28) together with that of the individual 
contributions given by (67) for one of them can be used to improve the form factors 
calculated in the hyperspherical formalism. In comparison with the most complete 8 
amplitude Faddeev calculation, the form factors, < S'|0|S' > and < S\0\MS >, obtained 
on the basis of Eqs. (61, 62), are found to be smaller by a factor 2-3 in the asymptotic 
domain around 50 (GeV/c)^. As to the ratio, ^^LL^^^ , an expectation can be obtained in 
the limit where one retains the dominant contribution (due to that part of the total wave 
function involving the terms obtained from Eq. (67) by exchanging the role of particles 1 
and 3 or 2 and 3): 

^<S\0\MS>, 1 k' ^ 

( <5|0|S> ''---72;rf<"-°-^^- <^^> 

It is found too small by a factor 1.5. Notice that the above factor in Eq. (69) is smaller 
than in Eq. (58) by a factor 4, showing once more the failure of the hyperspherical har- 
monic formalism to make accurate predictions in the present domain when it is restricted 
to the lower K values. As to the other form factors, < MS\0\MS > and < MA\0\MA >, 
a prediction made in the spirit of Eqs. (61, 62), but accounting for the minimal depen- 
dence of wave functions on k^ or K3, leads to a g~^° logg dependence of these quantities 
at high q. This implies that they should tend to zero with q more rapidly than the matrix 
elements, < 5*10 IS* > and < 5*10 IMS' >, while the full calculations rather indicate an 
opposite situation around q^ = 50 (GeV/c)^. The slow convergence to the asymptotic 
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behavior we do find in some cases, especially for the form factors, < MS\0\MS > and 
< MA|0|MA >, calculated in the hyperspherical harmonic formalism, invites to some 
caution in concluding from this discrepancy. 

In order to understand these results and their difference with the most complete ones, 
let us note that to get a q~^ asymptotic dependence of the form factors, one has to pick 
up terms proportional to p^ in the expansion of the product of the initial and final wave 
functions at small distances. This follows from dimensional arguments (the elementary 
volume involves six powers of p). To get such terms, one can combine p^ terms coming 
from the wave function in the initial (final) state with the zeroth order wave function (the 
wave function at the origin) in the final (initial) state. The contribution of these ones to 
the form factors are likely to be accounted for by Eqs. (61, 62). One can also combine 
terms of the first order in p in each state. These ones have obviously no contribution to 
the asymptotic expressions of the form factors given by Eqs. (61, 62). 

The failure of Eqs. (61, 62) to account for the asymptotic behavior can also be analyzed 
by considering the structure of the wave functions in momentum space. With a wave 
function as simple as (63), the dominance of the low momenta components in the matrix 
element (59), which is essential to get Eqs. (61, 62), necessarily applies to either wave 
function of the initial or the final state. With the wave function (67), which results from 
a first iteration over the previous one, but contains different factors, other possibilities are 
offered. The dominance of the contribution of low momentum components may involve 
one factor in the initial state wave function with respect to some integration variable and 
a factor in the final state wave function with respect to another integration variable. This 
is illustrated by the following contribution of the dominant part of a matrix element (Eq. 
(59)) involving the Faddeev amplitudes ■?/'i3,2 "^23,1- 



+ 



+ ( 



dKi2 dn^ 



(^2^ + ^2+^2)2 ((_i^3_^^^2)2+^2)2 



-W^ + J\^ + f '?12)2 + P^f (k?2 + {^Z + A/fg)' + 'P'Y. 



[^h + ^i + /9^)^ {^i + /32)^( {-Wz + Vfg) - f '^'12)^ + /3^)' 



.2'^12 ^ 2 '^SJ 



i4 + f'?;2 + vf^")^+/9T 



Ojt\j-\0 ClKjn 



■^12 



<2'4~/3 



«i + /.f +^2)2(4" +^T((-|('?3 



■ 1 ^" 
K 



2*^12 "T "'^S 



^/3^f"^2, ('_! 



^r'^.s) + ( o'^.i H o~'* 



v^; 



|g-)2 + ^2)2 



'q)-ff^i2y+^'r 



'Clrii n Clrio 



«;'2,4'~^ 



W 



( (i(/?'3 + 24' + ./|g-)2 + ^/2 + ^2)2 (^/2 + ^2)2 (^.2 + ^2). 
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^ ~ drio ClK/o 



(70) 



The two first terms on the r.h.s. of (70) are nothing but a contribution similar to that 
given in Eq. (60). They correspond to the processes represented in Fig. 9a. The third 
term involves the short range radial dependence of wave functions of both the initial and 
final states. It corresponds to the contribution of a diagram like that drawn in Fig. 9d. 

The integration in the above equation can be performed in the limit of large q, where 
the q~^ dependence factorizes out. The result is: 

where the two terms respectively correspond to the two first terms and the third one on 
the r.h.s. of Eq. (70). It is worthwhile to make a comparison with the contribution of a 
"diagonal term", ^^13,2 V^i3,2, or ^23,1 V'23,1, given by: 

- — . (72) 

It is seen that the contribution from a non-diagonal term, Eq. (71), is more efficient to 
contribute at high q, roughly a factor 50 larger. This is in relation with the expectation 
that the form factor at high q requires a momentum transfer from the quark which is 
stroken by the photon to the two other quarks. It is likely that a more complete calcu- 
lation will produce for the "diagonal term" an extra q~^ factor at least, which the above 
suppression factor partly accounts for (see Sect. 7.3 below). 
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Figure 9: Representation of diagrams contributing to the asymptotic form factor with the 
indication of various factors entering the gluon or quark propagators. These ones, which can be 
written as l/(n(|)^), are represented in the figure by 1/n, the factor (|)^ being factored out. 
Reading diagram a) from left to right for instance, one exphcitly gets: n = 1 for the first gluon 
propagator, n = 1 + 1 for the first quark propagator, n = 1 for the second gluon propagator 
and n = 4: + 1 + 1 for the second quark propagator. The weight accounting for undisplayed 
diagrams is also indicated: a factor 2 for similar ones where the role of the final and initial 
states is exchanged and, in most cases, a factor 2 for similar diagrams where the role of particles 
1 and 2 are exchanged. Factors accounting for the different spin-isospin components of the wave 
function are not incorporated here (rather see Table 2). 
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7.2 Further development 

From Eq. (59), and the wave functions provided by Eq. (67), one can now make a more 
complete calculation of the dominant term at high g, including some contribution due to 
the spin-spin force. The different form factors then read: 



<5|0|S> = ™^' 



aniqY 



n q^ 



<S\0\MS> 



729^6 (a m^ 



512, 1 1 2 + ^3, ,, ,,2 
^^ + — log -^) + 9 /t + <^ 

2 



^32, 1 1 , 2 + ^3, 15., ,,^ , 

— (^H log 1=) H ('« + /«^)2k'. 



729^31 a m„ 



i2 



< MS\0\MS > = '"7 "^^'7^ [0 + 3 ]{2k'^)\ 



TT q^ 



729v^ 



arrir 



\2 



< MA\0\MA >= - — ^-^--i^ [0-12 ]i2K'f, (73) 

TT g* 

The first term in the squared brackets corresponds to the improvement essentialy based 
on Eqs. (61, 62) discussed in Sect. 7.1. The second one corresponds to a further im- 
provement, which goes beyond these equations and was discussed at the end of the same 
section. The two terms have the same origin as the two first terms in the r.h.s. of Eq. 
(70) on the one hand and the third one on the other. They can be shown to have a close 
relation with the diagrams of Fig. 9 (see below). 

A slightly different but more complete expression, as far as the asymptotic behavior 
and the effect of the spin-spin force are concerned, can be obtained. We start from a 
zeroth order, Eq. (115), strongly peaked at zero momenta, and proceed to two iterations 
obtaining a wave function, Eq. (120), which is more appropriate than Eq. (67) as far 
as the short range correlation is concerned (no average over the angle), but is not 
normalizable. Then, the asymptotic behavior is given by : 

<S\0\S> = I^i^[7((. + <)^-20C + 9(«: + <)^] 
TT q° 

= (3.1C+17.9)(GeV/c)^g-^ 
= (19.8, 21.0) (GeV/c)%-^ 

TT g" 4 4 

= -(7.9C+13.2)(GeV/c)^g-^ 
= -(17.9, 21.1) (GeV/c)^g-^ 

<MS\0\MS> = 729v^(am,)^ ^ 3 ](2<)^ 

n q^ 

= 9.3(GeV/c)^g"^ 
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729^31 a m„ 



<MA\0\MA> = ^—^ — -^[0 - 12](2k') 



= -37.2(GeV/c)^g-l (74) 

The factor C accounts for extra non-perturbative effects, not included in the value of the 
wave function at the origin, a (see appendix C). Although there is no close relationship, 
differences between (73) and (74) can give some insight on the effect of the truncation of 
the Faddeev amplitude (remember that the wave function from (67) can be assimilated to 
the two amplitude Faddeev calculation). It is noticed that they only affect the first term 
of each bracket. The main difference is represented by the contribution of the term —2k,'^ 
to < 5*10 IS* >. This one could not occur in results presented in Eq. (73), which were 
obtained using the fully spatially symmetric component of the wave function, Eq. (63), as 
a zeroth order. The absence of difference for the second term in the bracket is due to the 
fact that, once the wave function at the origin, a, is determined, the calculation involves 
the first order terms of an expansion of the wave function around the origin. These 
terms, see Eqs. (113, 114), are completely determined by the - behavior of the potential 
at short distances. The numerical values have been calculated using a = 19.64 fm~ , 
mq = 337 MeV, k = 102.67 MeVfm and k,'^ = 1.66 n, the range given in the brackets being 
determined by the values C = 0.6 and C = 1.0. This range can be directly compared to 
the values that form factors shown in Fig. 8 are taking in the asymptotic domain. 

Table 2: Detail of the coefficients entering Eq. (74) with indication of the two interacting quark 
pairs and the diagram of Fig. 9 they correspond to. 
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V23VnO 
a 

Vi2{Vn + V23)0 



{V,3 + V23)VuO 

c 
sum 



V23OV,:, 8 

d 

(V^is + V2s)0Vu 1 

e 

sum 9 



16 _8 

3 3 



4 1 

3 3 



1 J_ 

3 12 



9 
"4 








-4 


4 


-12 


1 

4 


-1 





15 
4 


3 


-12 



33 



The different factors appearing in Eq. (74) can also be directly obtained, up to an 
overall factor, from considering the five types of diagrams shown in Fig. 9. The detail is 
given in Table 2. Three of them (a, b, c) involve two gluon exchanges in the initial or the 
final state wave functions. The two other ones (d, e) involve one-gluon exchange in both 
the initial and the final state wave functions. These two type contributions correspond to 
the first and second terms in the squared bracket appearing in Eqs. (74). They involve the 
wave function at the origin multiplied by its second derivative, ip^^\0) ip^^"^ (0), for the first 
one and the product of the derivatives at the origin, ^^^^' (0) ^^^^' (0) , for the second one. 
These expressions are those that are relevant for determining the asymptotic behavior of 
the form factor for the three-body case, while the expression, ^^^\0) %jj^^^ (0), is appropriate 
for the two-body case (incidentally, this expression works for the hyperspherical harmonic 
calculation hmited to the wave K = (see Eq. (53) together with Eqs. (37, 38)). 
We finally remark that, at high q, the combination of the various form factors, which 
corresponds to a scalar-isoscalar quark density, has the peculiarity to factorize into a 
simple form: 

< S\0\S > +i(< MS\0\MS > + < MA\0\MA >) = 



729^3 (am, ^2 



TT 



[7C + 9]{{K + K'^y-2K'J). (75) 



7.3 Remaining discrepancies 

The comparison of the Faddeev results in the asymptotic domain with the expectation 
given by the matrix element < 5'|0|S' > of Eq. (74) improves over that made previously 
in this section (7.1). On the other hand, the results of the eight amplitude Faddeev cal- 
culations are definitively closer to what is expected from (74) than the two amplitude 
Faddeev calculations. The remaining discrepancies may have their origin either in the 
truncation of the total Faddeev amplitude or in the fact that the asymptotic behavior 
has not been reached yet, unless both are related. At present, it is not clear whether the 
discrepancies should correspond to corrections of the order, -, as in the case of hyper- 
spherical harmonic results obtained previously, see Eqs. (55) and (56), or of the order, ^ 
(up to logg terms). Numerically, either one is possible in the higher range of momentum 
transfers we explored. 

From examining the general expression for the form factor, Eq. (59), or more partic- 
ularly the expression for the wave function obtained from Eq. (67), it is reasonable to 
consider that the theoretical corrections to the leading term should be of the order ^ (up 
to logq terms) in a complete calculation. A rigorous demonstration is difficult however. 
Any attempt to show it from an expansion of the form factor in powers of - indicates the 
existence of terms that linearly diverge in the limit g -^ oo, as it can be checked from the 
expression of a particular contribution, see third term in the r.h.s. of Eq. (70) (discarding 
the limitation k,'^, k,'^ r^ f3 ). In practice, they are cut off by a factor q, hence corrections 
of relative order -. It is possible that the total sum of these terms cancels out, leaving a 
correction of relative order ^, while in an incomplete calculation they would partly sur- 
vive. An example, though referred to the two-body case, is detailed in the appendix D. 
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It also provides some hint for making rigorous statements as to corrections of the relative 
order ^ to the leading order contributions. 

On the other hand, from QCD based calculations, the opinion is that corrections should 



be of the order ^ |31| , but this estimate incorporates the full quark propagator and, espe- 



cially, accounts for Z-type diagrams neglected here. One cannot infer from this calculation 
that the present one should evidence similar features. Notice that QCD based arguments 
hold up to logg terms. These ones partly originate from a relativistic treatment and 
should not be confused with those mentioned above, which are due to the three-body na- 
ture of the problem under consideration. The absence of convergence of some form factors, 
< MA\0\MA > in both the hyperspherical and Faddeev aproaches and < MS\0\MS > 
in the hyperspherical approach, makes the situation more confusing. As they are sensi- 
tive to the spin-spin part, one may expect them to depend on its short range behavior, 
which deviates from a Coulomb type potential we considered in our theoretical analysis. 
Assuming that the effect of the exponential factor in the force can be approximated by 
exp(-p/ro) with p = 3/q, part of the large dependence of the above form factors could 
be explained. No similar dependence is seen in the other form factors however. Further- 
more, one should not forget that we are dealing with a coupled channel problem, which 
can change the behavior of wave functions at short distances with respect to expectations 
based on perturbative arguments. In fact, this is a likely explanation for the slow con- 
vergence of the form factors, < M A\0\M A > and < MS\0\MS >, to their asymptotic 
behavior. 



7.4 Physical form factors 

Throughout this section and the previous one, we discussed particular contributions to the 
physical form factors with the idea to facilitate the analysis of their asymptotic behavior. 
Using Eqs. (74), the expressions of the physical form factors in the asymptotic regime 
are: 



729V3 {a rr ^^ 

8^ 



GU<f)- --vo(am. 



TT q 

{7C + 9)iK + K'^f - (9C + 15)(fi: + <) < - (14C + IS^^ 

= (-8.1C - 14.7) (GeV/c)« q-^ = -(19.6, 22.8) (GeV/c)^ g-^ 
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1 4 28 

i-C + 6){k + <)2 + (3C + 5){k + <) < + (-C - 20)<2 

= (-2.6C + 32.3) (GeV/c)^ q"^ = (30.7, 29.7) (GeV/c)® q'^ (76) 

The negative sign for G^^{(f), which imphes a change in sign for this form factor between 
and oo, is in agreement with numerical calculations presented in Fig. 6. This feature, 
which contradicts what is generally expected from QCD, has its origin in the unrealistic 
value of k'^ in the Bhaduri et al's quark-quark force (see next section). As expected from 
Eq. (76) above, the proton magnetic form factor is much smaller than the neutron one in 
the asymptotic regime. Again, this feature is strongly sensitive to the value of k,'^. 

Comparison of the different calculations of the form factors considered along the text 
with Eqs. (76) shows that the 8 amplitude Faddeev calculation is the most convergent to 
the asymptotic behavior. It cannot be said however that it has definitively been reached 
in the range 30 — 50 (GeV/c)^, confirming what can be inferred from examination of Fig. 



8 Making the calculations of the form factors more 
realistic 



In the previous sections, we mainly discussed the reliability of predictions with respect 
to technical aspects of the calculations regarding the description of the nucleon in terms 
of quarks. Any comparison with measurements was eluded. The reason for that is the 
existence of other effects to be accounted for. They are successively discussed in the next 
subsections. 



8.1 Quark form factors 



Y(q) Y(q) 




p,co 



a) b) 

Figure 10: Extra contributions to form factors that have to be included, due to a coupling to 
a quark-antiquark pair (a) or a vector meson (b). 
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In a constituent quark model like the one we used, it is quite legitimate to account for 
some electromagnetic form factor of the quarks. Contributions due to the coupling of the 
photon to a quark-antiquark pair, shown in Fig. 10a, or to a vector meson, uj or p, shown 
in Fig. 10b, cannot be generated by the non-relativistic constituent quark model and, 
consequently, their effect has to be added by hand @, Q. Processes represented in Fig. 



10b lead to a quite successful phenomenology at the nucleon (and meson) level, which is 
known as vector meson dominance. This can be transposed at the quark level without 
any change. The same statement partly holds for the contribution of the pion cloud to 
the proton and neutron form factors, which may be accounted for in some part in the 
p contribution. Some contribution may also come from the Darwin-Foldy term (for the 
charge density), but this one may be partly incorporated in the u or p contribution. 

For the proton charge squared radius, the latter can contribute an amount of 0.4 fm 
to be added to the matter squared radius, 0.25 fm^. The pion cloud, for that part not 
included in the p meson can contribute an extra 0.10 — 0.15 fm , making the proton 
and neutron charge squared radius close to the measured ones. At non-zero momentum 
transfers, the same mechanism tends to make the prediction for the form factors closer 
to measurements in the range, g^ = — 1 (GeV/c)^ [^. At higher momentum transfers, 
the hypothesis of the vector meson dominance doesn't lead to any improvement. On the 
contrary, assuming the full validity of this hypothesis, it adds an extra factor, g~^, to 
the discrepancy by a factor g~^ between the result of non-relativistic calculations and the 
QCD expectations. 



8.2 Improved description of the nucleon spectroscopy 

An other aspect concerns the force between quarks. The Bhaduri et al.'s force is known 
to miss the position of the Roper resonance by a large amount. To remedy this situation, 
one proposal based essentially on a phenomenological basis, has been to add a three- 
quark force |jTl|. This one has not led to any significant discrepancy in the low energy 



phenomenology of pionic baryon decays. On the contrary, it improves the situation in 
many cases [^, Q. From the expression of this force: 



2 



i^k^i^j^k rnoTij moVik 



which involves two Yukawa factors, it is expected that the form factor should have a q~^ 
asymptotic behavior, instead of q~^. This might contribute to make the form factors 
closer to the g~^ behavior expected from QCD. While the above three-body force may 
account in a phenomenological way for omitted relativistic effects, we nevertheless believe 
that the effect it provides has no relationship to the boost effects that could provide the 
appropriate behavior and should naturally appear in a full relativistic calculation (see 
Sect. 8.4 below). 

Calculations using the nucleon wave function obtained from solving the corresponding 
Schrodinger equation and the above three-body force have been performed in the hy- 
perspherical harmonic formalism with K = and K = 2. Due to a strong short range 
attraction (the force leads to a ^ term, similar to the centrifugal barrier), the solution 
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1 



at short distances, instead of 1 for the Bhaduri et al.'s 



behaves approximately hke 

potential [|14|. It results that the dominant contribution to the form factor, < S'|0|S' >, 
scales like q~^ at large momentum transfers, differing from the expectation by three pow- 
ers of q. The discrepancy has its source in the inability of the hyperspherical harmonic 
formalism with a limited number of terms to correctly reproduce the short range correla- 
tions between quarks produced by the three-body force given by Eq. (77). The problem 
is similar to that encountered for the Coulomb part of the force in the same approach, 
with the difference that it is more severe here. The calculated form factor is shown in 
Fig. 11 (dash-dotted line). 




100 



Figure 11: Form factors calculated with three-body forces and without. Those presented here 
correspond to the < S'lOIS' > transition matrix element. They have been multiplied by a factor 
q^ to evidence the convergence to the asymptotic value expected for two of them (potential of 
Bhaduri et al. and potential including a gaussian three-body force, Eq. (78), continuous and 
dashed lines respectively). The other curves represent results for the three-body force given by 
Eq. (77) (dash-dotted line) and for the Bhaduri's model in absence of spin-spin force (dotted 
line). 

Due to the proximity of some collapse in calculating the nucleon wave function with 
the above mentioned three-quark force, calculations were also performed with a Gaussian 
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rp^ 



Vjfj = Voexp{-Y.'-i^)- (7J 



type force: 

v}?j = Voexpi-J2- 

This force is much less singular at short distances than V^y and tends to zero more quickly 
for high momentum transfers, like exp(— ^(/t^g + ^i))- The resulting nucleon spectroscopy 



is equally good [Q. The form factors so obtained do not differ much from the previous 
ones and, like them, overshoot those obtained with the Bhaduri et al.'s force by a large 
factor in the range around 50 (GeV/c)^ (see Fig. 11, dashed line). Contrary to the other 
calculation, the present one roughly evidences a behavior g~^. This result is expected 
in the hyperspherical harmonic formalism with a restriction to -ft" = and K = 2 when 
the Coulomb part of the potential takes over the gaussian part. Surprisingly, it does not 
seem to be strongly affected by the large strength of the gaussian three-body force. While 
the asymptotic power law is identical to that obtained with the Bhaduri et al.'s force 
(continuous line in Fig. 11), the overall size is changed by three orders of magnitude. 

The comparison of the various results obtained from forces without and with three- 
body forces tends therefore to show that extra contributions to the description of the 
force between quarks, especially at short distances, can considerably modify the onset 
of the asymptotic behavior of the form factors, making it to occur at higher momentum 
transfers or with a different pattern. Thus, the form factors in the model involving the 
three-body force given by Eq. (77) reach their asymptotic behavior from below. As to 
the form factor calculated with the other three-body force, Eq. (78), the apparent onset 
of the asymptotic behavior that is seen in Fig. 11 is misleading. The product of the 
form factor multiplied by g^ represented in the figure passes through a maximum and 
then should decrease to reach a value which is of the order of 0.710^ (GeV/c)^. Except 
for a difference in scale, the pattern exhibited by this result is quite similar to that one 
obtained for the Bhaduri et al.'s force without spin-spin force (dotted line in the same 
figure). 

We would also like to notice that the difference by three orders of magnitude for the 
form factors at high q produced when three-body forces are included is not a distinctive 
feature. A similar effect |^^ has been obtained with a two-body quark force resulting from 



meson exchange [27|. The common feature has perhaps some relationship to the small 



mean squared radius of the nucleon that these models predict, of the order of 0.10 fm 



8.3 Description of the spin-spin force 

It is known that the spin-spin part of the quark-quark force we used, essentially fitted on 
the difference of the vr- and p-meson masses, is too large with respect to a genuine gluon 



exchange force by a factor 6 [33]. Correcting for this factor would decrease the contribution 
of the mixed symmetry states to the form factors and, indirectly, would accelerate the 
convergence to their asymptotic value since these particular contributions have a slower 
convergence. It would also probably remove some of the sign changes that are observed 
for some calculated form factors, especially G^^{(f) and G\^{(f). These changes are not 
expected to occur and, furthermore, are not seen in the measurements. 
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The correction by a factor 6 is not the main problem however. Apart for the fact that 
it should be compensated by something else, the radial part of the force, which is given 
by a (5 function in the non-relativistic limit, has been replaced by a less singular function, 
namely a Yukawa type potential. The advantage of this one is that it can be incorporated 
into a Schrodinger equation while avoiding a collapse of the solution in some cases ||3^ . 



It is noticed that the spin independent part of the force also contains a S type force, but 
this one is generally neglected. 

In a perturbative calculation of the asymptotic form factors, there would not be any 
difficulty to use a 6 type force. The propagator of the gluons in Figs. 2 would be replaced 
by a constant, removing for each of them a g~^ factor in predicting the behavior of the 
asymptotic form factor. The result so obtained should be closer to the asymptotic QCD 
prediction. This feature however supposes approximations that are not supported by a 
more rigorous derivation of the force. The 6 type force is expected to be made less singular 
by the presence of relativistic normalization factors, ^ ^. The resulting non- locality is 
difficult to account for, and in absence of a better solution, employing a Yukawa potential 
instead of a 5 type avoids to make unrealistic predictions. 



8.4 Relativistic corrections 

A complete relativistic calculation of the form factor is rather tedious. Some recipe to 
account for part of these effects has been used in the past. One of them consists in 
changing the argument of the form factor as follows |3^, ^, ^ : 

F-{q') = F-{-^). (79) 

This expression, which has its origin in a naive estimate of the effect due to the Lorentz 
contraction, leads to a constant at large q, in contradiction with the QCD expectation. 



It turns out to be valid only at small q ^^ |4^. A modified expression, motivated by the 
QCD result, has also been used ETI, 



The discussion presented in Sect. 2 indicates that the non-relativistic calculation misses 
an important factor due to a boost effect which mainly affects particles interacting by 
exchanging spin 1 bosons. This factor may be roughly of the order of: 



' 2mq ' 



^0) 



where Eq is an energy corresponding to an average momentum, g, carried by the quarks 
in the Breit frame, see Fig. 3. At large q, this energy varies like a fraction of q and the 
above factor, Eq. (80), provides the q^ factor discrepancy between the non-relativistic 
and relativistic calculations, while allowing one to recover the non-relativistic prediction 
in the limit of small momentum transfers. There is no point to insist on the approximate 
character of this factor. It results from an analysis of what a non-relativistic calculation 
partly misses and the determination of the coefficient of the g^ factor in Eq. (80) is likely 
to be uncertain. Furthermore, corrections of the same order are expected from Z-type 
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diagrams, which involve the negative energy part of the quark propagator. In absence of 
a more refined analysis, we will simply mention an estimate made by incorporating in our 
calculation the factor given by Eq. (80) together with a value of q between | (momentum 
of initial and final quarks in the Breit frame) and | (momentum of quarks in the interme- 

4 

diate state of Fig. 9.d). This amounts to a factor ^^^^ ^4 at high q. Disregarding the effect 

of the mixed symmetry component due to its uncertain character (see above), we consid- 
ered the contribution of the matrix element < S'lOjS' > to the quantities q^ G^j^cf) / ^p 
and q'^ G^ji^cf) / iin- Results so obtained for the Bhaduri et al.'s quark interaction model 
have the size that experiment suggests, 0.4 (GeV/c)"^ [^, but are off by three orders of 
magnitude for the other models incorporating three-quark forces. 



9 Conclusion 

We have studied the nucleon form factors using a non-relativistic constituent quark model. 
The emphasis has been put more on general features pertinent to their calculation than 
on a fine description of observables. A special attention has been given to the asymptotic 
behavior of the form factors at high momentum transfer. 

Contrary to what is sometimes mentioned, it is found that the non-relativistic quark 
model does lead to power law form factors at high q. This power law is related to the force 
between quarks in momentum space and thus can be predicted if it is well defined. It is 
our belief that this approach, possibly corrected for exchange currents and relativity, is 
completely equivalent to the one used in the QGD asymptotic regime where uncorrelated 
wave functions are used together with an effective operator accounting for three-quarks 
absorbing a virtual photon while exchanging two gluons. In the present approach however, 
the calculation of the form factors is performed in a unique model of the nucleon and tends 
to cover the full range of momentum transfers from to 00. 

The exact power law behavior, which is expected to be q~^ in the non-relativistic frame- 
work used in the present study, has been found to depend on the approximations made 
in treating the three-body system. While the hyperspherical harmonic formalism with a 
restriction to the lowest values of the grand orbital, K, is generally a good approximation 
for determining the binding energy, it has been found to systematically lead here to the 
wrong power law, q'"^ instead of q~^. To remedy this situation, the full set of K values 
is required. This is expected from the behavior of the wave function of two quarks at a 
short distance, on which it depends linearly. As far as we can see, the Faddeev approach 
is unsensitive to this difficulty, but we did not get the certitude that the magnitude of 
the asymptotic form factor is independent on the truncation of the Faddeev amplitude 
although it could be. In any case, this dependence seems to be small. At low momentum 
transfers, the form factors calculated with the 2 amplitude Faddeev approach or with 
the hyperspherical formalim are in good agreement with the more complete 8 amplitude 
Faddeev calculation up to g^ = 2 (GeV/c)^ for the former and g^ = 3 (GeV/c)^ for the 
latter. 

The convergence of the form factors to their asymptotic value is another feature 
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we looked at carefully. By only considering form factors in the range around q^ = 
50 (GeV/c)^, it has often been difficult to guess their power law behavior. This has led us 
to examine the wave functions at short distances and, from them, to determine the asymp- 
totic behavior. In a few cases, we had to conclude that form factors at q^ = 50 (GeV/c)^ 
were off their asymptotic value by a factor 2-3. This occured especially for form fac- 
tors involving the transition between mixed symmetry components, MS <-h> MS and 
MA ^-y MA and is due in these cases to sizeable next to leading order contributions, 
which seem to be of the relative order q^^ but may also be of relative order g~^ logg as 
generally expected theoretically. 

Most often, the behavior of the form factor at high q is illustrated by a process where a 
quark, which has been stroken by a virtual photon, shares the momentum it received with 
the other two quarks by the successive exchange of two gluons. In practice, calculations 
involve several processes with different time orderings of the photon absorption and gluons 
exchanges. In the present work, the dominant contribution rather comes from a process 
where the same quark that is stroken by the photon has previously and subsequently 
exchanged a gluon with each of the other quarks. It is likely that this should also hold in 
QCD, but this has to be checked. In any case, we believe it is worthwhile to emphasize the 
point as some arguments and calculations have been developped without regard to this 
type of process. Its contribution is essentially determined by the square of the derivative 
of the wave function at the origin. This differs from the two-body case or the simplest 
calculation with the harmonics hyperspherical formalism, which only involve the first 
derivative. 

While we have not been able to fully answer some of the questions raised in the introduc- 
tion, we showed that some caution is required in dealing with the asymptotic behavior of 
form factors. Their onset may not occur as quickly as usually expected. Most important, 
it supposes an accurate determination of the three-quark wave function at short distances. 
Present calculations provide enough evidence that the high momentum transfer behavior 
in a non-relativistic approach is given by q~^. We believe that the difference with the 
power law, g~^, expected from QCD is most likely due to relativistic (boost) effects that 
have not been considered in the present work. Consequently, it is inappropriate to use the 
g"'^ behavior to get constraints on the non-relativistic three-quark wave function. Con- 
trary to a belief that probably originates from a non-relativistic approach or from studies 
with spinless particles, a g~^ behavior of the form factor does not imply a similar behavior 
of the wave function at high momenta. Along the same lines, it is inapproriate to use a 
truncated hyperspherical harmonic approach to analyze form factors at high momentum 
transfers. In either case ||2^, ^ |13|, |3^ , a bias is systematically introduced, which prevents 
one to make relevant conclusions. 

For the future, it would be worthwhile to study more carefully the next to leading order 
corrections to the form factors, both mathematically and numerically. The questions to 
be answered are the relative order of these corrections, -, or ^ (up to logg terms), and 
the sensitivity to the truncation generally performed in solving the Faddeev equation. 
Obviously, a full relativistic calculation would be desirable. It is likely that the questions 
raised here will appear there too. We hope that the answers we got will be of some 
relevance for this more difficult problem. 
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A Wave functions in the hyperspherical harmonic 
formalism 

A.l Jabobi coordinates 

For a wide range of physical potentials (mainly two-body potentials depending on the 
relative coordinates), the three-body problem is conveniently treated by using Jacobi 
coordinates. For the equal mass {m) case, the Jacobi coordinates can be defined as 
corresponding to the vectors: 



ri-r2 

Pl2 = 



V2 ' 



P3 



^S{f,-^^). (81) 



Altogether with the center of mass (i?) coordinate, they determine completely the position 
of the system. 

The kinetic energy operator can be separated as 

T = y^ = — + ?j^ + l£^^T,ra+T,nU (82) 

2m 6m 2m 2m 

t 

containing the center of mass kinetic energy, ^ = Tcm.-, plus the internal kinetic energy. 
Tint, in terms of the Jacobi momenta. Then for potentials depending only on the Jacobi 
coordinates, the center of mass motion can be factorized out. 



A. 2 Hyperspherical coordinates and harmonics 

From pi2 and p3, one can define hyperspherical coordinates p G [0, oo), G [0, 7r/2] in the 
form: 

Pi2 = p sin 0, 

P3 = pcos0, (83) 

the differential volume element d^pi2 d^Ps being then written as: 

p^ dp dVt = p^ dp d(j) sin^ cos^ dpi2 dp^. (84) 
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The internal kinetic energy is then expressed in terms of p and Vt = {612, V'12, ^3, V^s, 4>}, 
(where {6, if) are the spherical angles of pu and ps) in the form: 



T^, 



int 



I ( d^ 5d_ K\ny 



(85) 



2m \ dp^ p dp p^ 

where K'^{Vt) is an angular operator whose eigenfunctions are called the hyperspherical 
harmonics (HH), Y\^k]{^)'- 

K\Q) V](fi) = -K{K + 4) V](fi), (86) 

K, the so called grand orbital number, defines the parity of the HH as {—)^. One can 

choose appropriate combinations of hyperspherical harmonics with definite values of the 

total orbital angular momentum L{L = I12 + /s) and its third projection M^ and with 

definite spatial symmetry: 

\ 

[(K,sym)] 



rit-'^'^J. m. 



In terms of these, the wave function of a system of total momentum P = and spin 
and parity, J^, reads 

\^>{p,n)>=<pn\m>= Y. i'KAp)^lKX:;\{msyrn>\^^^Y.^,{p)Y,{n), (87) 

sym 

where \sym > stands for a spin-isospin wave function of definite symmetry and the square 
bracket implies the coupling to total angular momentum J with the required overall 
symmetry. The subindex j of ipj{p) and Yj[Q) corresponds to a simplified notation, j 
comprising K, angular momentum and spin-isospin variables. 



The state |\E'(p, ^2) > satisfies the Schrodinger equation: 

+ -7^ + TT^ ] +V + 3m 



2m \ dp^ p dp 



E 



p- 



|^(p,f]) >=0. 



By substituting |\l'(p, Vt) > by its expansion, Eq. (87), and projecting it, one gets a system 
of coupled equations (its number depending on the number of Yj{Q) terms kept in the 
expansion) : 



2m \ p^/^ dp' 



;P 



5/2 



{K + 2)^-1/4'' 



+ 3m-E 



i^jip) 



+ J2Vjf{p)i'f{p) = o, 



^9) 



where the matrix elements of the potential are given by 

Vjj,{p) = f dn Y*{n) v{p,n,s,T) Yj{n), 



(90) 

In this equation, s, r indicate the possible dependences of the potential on spin and isospin. 
The normalization condition is: 



5:/rfpp^iv^,(p)p 



1. 



(91) 
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A. 3 Details about wave functions 

For the moment, we are going to truncate the expansion (87) by considering only K = 0,2 
terms. Furthermore, if we consider a restricted potential as the one given by (17) but 
without spin-spin interaction, the ground state does not contain K = 2 component. Then 
we can write: 

\^ip,n)>,+=^Mp)\S>, (92) 



Mp) = 0. (93) 



TTVvr 
satisfying 

1 , 1 rf2 ,,„ 15, AV2k 16^2 p 3D 



2mg p^/2 dp^ 4p2 T\ p hix d^ 2 



The terms — ^^- and ^4^^ represent respectively the monopole components of the 
Coulomb and linear parts of the quark-quark interaction: 

-l^!!._3/dnf ' )^ / . (94) 

TT p J yn^-n j 2v2psin0 



16x/2p , ( ,^( 1 \ pV2 



sm 



3 \d^\^\ ^^^— ^. (95) 



57r a^ J \'n'\/TT I 2a 



The relative distance ri2 appearing in the potential (17) has been replaced by its expression 
in terms of the hyperspherical coordinates, \/2p sin0, and the factor 3 accounts for the 
interaction between the three pairs 12, 23 and 13 which contribute equally to (94, 95). 
The solution of Eq. (93) at short distance can be expanded in terms of the powers of p. 
It reads: 



^i(p)p^o = a'i 



8^2 

— — Kirig p 

oil 



1 / , 3D, ,8^2 ,n\ 2 

~i2 r"^«*^^ ^ ^"^'? ^^^~ ^(^^r'*"^^'' J ^ 

1 /32^m„ Q8V2 ,^ ^ 3D, 2 25, 8v^ ,3^ 3 

21 V 57r a2 157r ^ <? ^ 3 ^ " 12^ hn 'i> j >' ^ 

where a'^ is a constant (the value of the wave function at the origin). 



(96) 



The coefficient of the term linear in p is determined by the Coulomb part of the force, 
while the energy of the system appears in the second term in p^ of the expansion. The 
intensity of the confining potential affects the third term in p^ of the expansion. It however 
implicitly appears in the energy E. 

A particular case of Eq. (96) is that one where the confining potential is neglected. 
Equation (93) then reduces to that of a Coulombian type problem. The solution is given 
by: 

^(p) = ae-^P, (97) 
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8^2 



64 2 



with P = ^nrug and E - Srug + ^ = -£^ - 

From the normahzation condition of the radial wave function given by: 



{91 



one gets a 



ml 
VT20' 



In the case where the spin-spin interaction is turned on, one has to solve a set of coupled 
equations. For our purpose here, it is enough to keep the most singular part of the spin- 
spin force (^ (Tj.(jj) (the complete expression has been used for the numerical calculations 
detailed in the main text): 



1 d^ 



2mn 
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(99) 
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-— + Sniq - E 



Up) = ^<-Mp), (100) 

on p 



where ipsip) = ipK=2,L=o{p)- In the limit p — > 0, it is easy to show that to remove the \ 
singularity on the left hand side of the second equation, ip^ip) has to go to zero as : 
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(101) 



Similarly, other components in the limit p ^ may be determined. An expression of 
the wave function beyond i^ = 0, 2, i.e. including K = A and K = Q terms, and limited 
to those terms linears in the variable p is the following: 
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(102) 



B Details about calculations of the form factors 



The complete expressions of the form factors calculated from the analytic expressions of 
the wave functions given by Eqs. (37) and (41) are as follows: 



< S\0\S >= A05Vq 



al/3i 



(g2 + 6/32)7/2' 



(103) 
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2 

< 5|0|M^ >= -2835v^ ^iS r, (104) 



a 



2 



/93V^ 



< M^|0|M^ >= 7290^ [QA - ^^, ;y^2)n/2 (105) 

175g^° + 574g^/?| + 924g^;3^ + 792g^;3^ + 352g2/3^ + Upf 



< MA\0\MA >= -2430^ ^64 - ^^, ^;^^,,,, (106) 



168g''' + 588g*/?| + 924g'^/5| + 792g^/3| + 352g"/5| + 64/^3^ 



where /^a = v^/^s- While close expressions of form factors can be obtained in a few cases, 
as above, quite generally and especially for more realistic wave functions incorporating 
terms depending linearly on the variables, r^, ri3 and r23, this is not possible. 

One may however be interested to get some prediction for the asymptotic form factor. 
Starting from the expression of the matrix element, < X|0|y >, given by Eq. (49), and 
after integration over the orientations of the pu and ps variables, which can be performed 
in most cases, one is left with a quantity of the following form: 

/(g) = — \ — d(f)sm'^(f)cos^(f) rfpp'' sin(W- gp cos 0) iJ(p, 0), (107) 

TV \j 2 q J V 3 

where H{p,(f)) involves the nucleon radial wave functions. The overall factor, — , arises 
from a factor (47r)^ due to the integration over the angles relative to the vectors, pi2 and 
P3, and another one, (^^A^)^, appearing in the normalization of the wave function, Eq. 
(102). Assuming that H{p,(f)) can be expanded at small p as: 

iJ(p,0) = Fo + pFi(0) + ... +p"F„(0) + .., (108) 

it is clear, from dimensional arguments, that the p" term will produce a contribution to 
the form factor proportional to g~(^+"), besides S{q) functions or its derivatives which 
are irrelevant at high q. To get the corresponding coefficient, a„, one has to perform the 
integration over the p and variables: 



16 6 + n/3 

V V2 



lim^^o / d(f) sin^0 cos^0 dx x"^"*"" sin(x cos 0) F„(0) e ^^. (109) 



As in Eq. (13), the factor e ^^ allows one to get rid of the part involving the S{q) function 
or its derivatives. After integration over x, Eq. (109) becomes: 



1^ %"/|lim,_o(-^)^+" /rf0sinVcos^0-|^ 
TT V 2 de J cos^0 + e^ 



As a check, one can show from expression (37) that 



i^n(0) = «?m^. (in; 
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Putting it in (110) and using the relation: 

, sin^(/) cos^(/) TT /I 



cos20 + e2 2 V2 



e 



-evOT^ , (112) 



the different terms of the expansion of the form factor (103) in terms of the powers of -, 
-sTTT {n = 0, 1, ..), are recovered. 



C Short distance behavior of the wave function in the 
Faddeev formahsm and asymptotic form factors 

C.l Short distance wave functions 

In the Faddeev approach, and for the considered interaction given by (17), the short range 
behavior of the wave function (28) can be obtained from Eq. (26). So, for the amphtudes 

^12!3(^l5^2,^3): 

3 

^12,3(^,^2, fg) ocl - - (k + 30 rUg ri2 + ..., 

3 

^^2,3(^15 ^2, fs) oc 1 - - (k - O rUg ri2 + .... , (113) 

and for the full wave function and up to an undetermined factor: 

^F,o(n,r2,r3) oc M - -(/t + K^) mg {ru + ri3 + r23) + ■■■) \S > 
-i<m, (((2ri2 - ri3 - rsa) + ..)\MS> +V3 ((r23 - ^3) + ..)\MA>). (114) 

It is worth to notice that from this wave function one can recover the first terms of the 
expansion, Eq. (102), by making the projection on appropriate hyperspherical harmonics. 
It is also noticed that this wave function, after averaging over all directions of fi2, has a 
term linear in its modulus, ri2, or equivalently in sin0, while the truncated wave function 
given by Eq. (102) depends on the square of sin0. As this linear term in ri2 corresponds 
to a first order in the interaction, the wave function in momentum space given by Eq. 
(67) should have some track of it. In fact by taking the Fourier transform of Eq. (67) it 
is possible to recover (113). The proportionality factor is then found to be k:;^/^- 

An expression for the momentum space wave function reproducing the short range 
behavior of the r-space wave function (114), may by useful. Starting from the following 
zeroth order completely symmetrical wave function (see permutation relations at the end 
of this appendix) , 

^S3 + *Sl + ^^2 = ^^(2^)' ('^('^3) S{k,2) 

+'^(-2'*3 + ^'^12) '^(2'*i2 + -Y^^' 
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+^(-2'«3-^/ti2) d(-/ti2- — Ks) I |5>, (115) 

one can iteratively obtain a first order wave function from tlie momentum space Faddeev 
equation. For one of the amplitudes, ^u^si^-u, ^,3) for instance: 

(^gy^, 4) + (1 -- 3) + (2 ^ 3)) , (116) 

where the potential, here given for definiteness for its Coulomb part and in accordance 
with Eq. (17), is defined as 

v/(/«i2,/ti2) = 



'?12 - /?12)^ 



Then 



^Ss + ^2!! + ^§2 = ^^m,C(i) (J3 [{k + k'^)\S > +2k'JMS >] 

+ Ji [(/. + <)|5 > +2<(-l|MS > +^\MA >)] 

+ J2 [(/. + <)|5 > +2<(-i|M5 > -^\MA >)]] , (117) 



where: 



h 



{2nf 6{kA 
'' («:?2 + /9«^)'^?2' 
(27r)3<5(-i/?3 + f/5i2) 



((i/5i2 + f /^3)2 + /3«')(i'^i2 + f %)2' 
(27r)3 ^(-|/?3 - fKu) 

((|/?12 - f /^S)^ + /?«')(i/5l2 - f /^S)^' 



:ii8i 



and C'^^^ is a normalization constant. The quantities Ji and I2 are obtained from Is by 
performing a circular permutation and re-expressing the variables in terms of k,i2 and ^3. 

The -A- behavior at high values of /ti2 is essential to reproduce the linear dependence in 

'^12 

fi2 in Eq. (114). 

From Eq. (116), multiplying both members by the factor K,'l2, taking the limit /ti2 -^ 00 
and integrating over the variable /?3, one gets for the Coulomb type potential: 



lim /ti2 / C?/?3^12,3(/«12,/t3) 

12-+00 ^^ J 



= ^ra,{K + <) J ^^^ (^i2,3(«:'i2, 4) + (1 - 3) + (2 ^ 3)) . (119) 

Up to a common factor, the l.h.s. represents the coefficient of the ri2 term (see the 
first line on the r.h.s. of (114)) while the r.h.s. represents the configuration space wave 
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function at the origin, \l/(0). Now, by using the wave function given by (117) and under 
the requirement of recovering the coefficient of the hnear term in r in(114) one gets 
P^^^ = -j=mq{K + k'^). In this way, the corresponding wave function can be used in 
calculations where this property may be required as in the calculation of the form factor 
in the asymptotic regime. This value also ensures that the wave function at the origin is 
unchanged at this order, with the result C-^-* = 1. 

It is also possible to iterate over the solution given by Eq. (117) ((114) in configuration 
space). The amplitude, \E'i2,3, so obtained, reads: 



12,3 
J \2/ J , T I T \ , /o. ' \2 



1 . 1 



[k + <)^( J3 + Ji + J2) + (2<)'( J3 - i^Ji - 2-^2)] \S > 
+ (/€ + 0(2<)[(J3 + Ji + J2) + (J3 - ^Ji - \j2)] \MS > 
+^('^ + 0(20(^1- ^2) |MA> 

+%2<)Vi - J2) \A >\ , (120) 



where J,- is related to /,■ as follows: 



with the result: 



^ 1 /■ dK[^ dn'^ 6{k3 - 4) , 

■'' kI, + 4 + /3(2)^ J (27r)3 (^12 - A2)' ' ' 



1 1 ^ ./3«, (27r)35(/53) 



[\1\\ 



J 



«?2+/5(')' «i2/9W ^«:i2' 47r/5« ' 
1 1 ^2 1 



^"-3; VO V3"'3^P- ^3' 
1 123 1 

^^ - ki, + ki+ p(^r (,5i2 - ^/53)^ ^73^ (M+^(^nFI' ^^^^^ 

For simplicity, we omitted the dependence of the functions, Jj on the arguments, ^[2 ^^^ 
Kn, what is reminded by a " ' " at I. 



12 

3' 



In the limit where P^^' can be considered as an infinitesimally small quantity, the last 
factors in Ji and J2 tend to the function S{ks), making the quantities Ji, J2 and J3 
equal to each other. Equation (120) simplifies and identifies to that part of Eq. (117) 
corresponding to the same Faddeev amplitude, \l'i2,3, proportional to I^. On the other 
hand, the presence of the S{k3) in Eq. (122) has no physical foundation. It results from 
the iteractive character of the calculation and should disappear in a complete one, as 
partly realized when going from the ffist iteration, Eq. (117), to the second one, Eq. 
(120). 

The parameter, P^'^\ may be determined by requiring that Eq. (119) is satisfied when 
the expression of the wave function, Eq. (120), is used. This can be done once P^^^ has 
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been determined from the previous iteration {P^^^ = -j^niq^K + k!^)). Disregarding this 

value, one may also assume /9*^^^ = Z?'-^-' and determine the corresponding common value 
of these parameters by fulfilling Eq. (119). The value of C*^^-* is determined so that the 
configuration space wave function at the origin, ^1^(0), is equal to that one in the previous 
iteration or its value from a full calculation: 

"^(^^ = / ^1^ (^12,3(42, 4) + (1 - 3) + (2 ^ 3)) . (123) 

For a pure Coulombian problem and depending on the approach, values C^^^ = 1.0 
and C*^^^ = 0.66, corresponding respectively to (3^'^^ = 1.105 Kirig and Z?'-^-' = 1.403 KTTig, 
are obtained, providing some uncertainty range. Ultimately, after an infinite series of 
iterations, the value of the parameter corresponding to P^"^^ should converge to that related 
to the binding energy. This one is expected to be equal to /5 = ^^nrriq in the hypercentral 

approximation {K = 0), Eq. (97), and /? = 1.016 ^^nniq = 0.732 nniq in a more complete 
calculation ET 



C.2 Asymptotic form factors 

Concerning the asymptotic behavior of the form factors, it is interesting to apply Eqs. 
(107-109) to the first order term of the wave function given by (114) which involves terms 
ri2 and ri3 + r23. After integration over the various angles and up to some factor, one is 
left with the following expression of -Pi(0): 



Fi(0) = sin0 (forri2), 

(forri3 + r23). (124) 



|sin0 + \/3cos0P — |sin0 — y^cos^P 



6v^sin0 COS0 
The first term in Eq. (124) can be dealt with easily. Using the relation: 

, sin^0 cos^0 . , 2 /, 9x /TT , /.„^x 

it is seen that, in absence of e' terms in the expansion, there is no contribution to the 
coefficient of the g~^ term, Oi. The second term in Eq. (124) can be written as a term 
similar to the previous one: 

/",, sin^cosV 1 vr 

/ ^'^ 2^ I 2 — — = l-nTT-arctge), (126) 

J cos"^0 + e^ sm0 2 

plus another one integrated over the range, | < cos0 < 1. As this one safely converges 
when e — i> 0, it is immediate that its expansion in terms of powers of e will exhibit only 
even powers and, therefore, will not contribute to ai. Thus, the linear term in p appearing 
in the expansion of the correct wave function around the origin does not contribute to the 
form factor at the order g~^, as expected. This is not a trivial result however. A linear 
dependence of the wave function on the variable p3, for instance, which has no physical 
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ground, but turns out to be undistinguishable from a contribution of a pi2 dependence 
in the approximated calculation retaining the K = wave, leads to a g^'^ asymptotic 
behavior of the form factor. This indicates that the pa dependence of the wave function 
has to be carefully determined, otherwise it would lead to a bias in calculating form factors 
in the asymptotic domain. 

Expression (120), which involves the effect of two gluon exchanges without any restric- 
tion, in contrast to Eq. (67), can be used to determine the asymptotic behavior of form 
factors. When it is averaged over the various angles, including 0, it allows one to recover 
the high momentum behavior of the wave function (63) obtained in the hyperspherical 
formalism, with the same front factor, independently of the value taken by C^"^^ . 

Throughout this paper, and especially to get the expressions given in this appendix, we 
employed the following permutation relations: 
exchange of particles 1 and 3: 

1^ VS^ -. 1^ V3^ 

P3 -^ -^P3 + ^Pl2, Pl2 -^ ^Pl2 + ^P3, 

\MS>^--\MS>+—\MA>, \MA>^-\MA>+—\MS>, (127) 

exchange of particles 2 and 3: 

1^ V3^ -. 1^ 73^ 

P3 -^ ~2^^ 2~^^^' ^^^ ~^ 2^^^ 2~^^' 

1^ V3^ ^ 1^ V3^ 

\MS>^--\MS> - — \MA>, \MA>^-\MA> - — \MS> . (128) 

I 2' 2 ' ' ' 2' 2 ' ^ ^ 



D Details about the q ^ asymptotic behavior of the 
form factors and terms beyond 

We first consider in this appendix the effect of a truncation of the total wave function 
calculated in the Faddeev approach, which has been limited to the ten partial waves with 
the lowest values of the angular momentum of the pair of particles 1 and 2: / = 0, 1, 2, 3 
and 4 for each spin state of the same pair, and 1. We especially look at the contribution 
of the term which arises from the wave function given by Eq. (114) and involves the 
quantity ri3r23. Its contribution to form factors, which is the dominant one, involves / 
values ranging from to oo, while the contribution of the other term ri2(ri3 +r23), which 
is less important, only implies / = 0. The first contribution may thus be sensitive to 
some truncation, while the second one cannot. Let's remind that the truncation under 
consideration is different from that on the Faddeev amplitude, for which two cases with 
2 and 8 partial waves have been considered in this work. 
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Using an approach similar to the one presented in the appendix B, we calculated the 
contributions to the coefficient of the q~^ factor due to different waves. We separated 
those for even and odd values of /, corresponding respectively to the matrix elements 
< S\0\S >,< S\0\MS > and < MS\0\MS > on the one hand, and < MA\0\MA > 
on the other. Assuming that the total sum is normalized to one, we found that the 
contribution for a given / is given by: 



(2/ - 5) (2/ - 3) (2/ - 1)(2/ + 3) (2/ + 5) (2/ + 7) ' 

The sum thus reads for even values of /: 

4 100 4 

1 = 0.5714(= -) + 0.4329(= — ) - 0.0040(= ^^) - ... (130) 

and for odd values: 

20 20 

l = 0.952(=-) + 0.046(=— ) + .... (131) 

The above result shows that no serious discrepancy is introduced by the truncation of 
the total wave function we made and that the origin of the discrepancy of the calculated 
form factors and their expectations in the asymptotic regime has to be looked for else- 
where. On the contrary, calculations using phenomenological wave functions limited to 
/ = would in any case miss a sizeable contribution. 

We mentioned in the text the possibility that the correction to the dominant q^^ term 
be of the relative order g~^, with the result of delaying the convergence of the form factors 
to their asymptotic value. We here give details about an example dealing with a two-body 
system. This is not directly relevant to the three-body case of interest in this work, but 
we nevertheless believe it could cast some light on this one. 

The expression of the form factor for a system of two equal mass particles is given by: 

F{t) = \^^^^{k)^(k^-^q), (132) 

where v^(A;) represents the wave function of the system under consideration. For the lowest 
state of the hydrogenic system discussed in Sect. 3, •^(k) = V^np^f;^- Hence: 

^ f dk GAttk^ k^ 



(27r)3 (F + n^)2^(k + ig-)2 + ^2)2 {±q2 + ^2y 
In the limit of large q"^, the form factor reads: 



i^(g^) = ^(1-32^ + ...), (134) 



'2 



evidencing a g correction to the leading order term. It is interesting to compare this 
result with that of a slightly different expression: 

^^ ^ " i ^2^ (A;2 + fi:2)2(A;2 + ig2 + /,2)2' (135) 
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corresponding to the first term of an expansion of the integrand in Eq. (133) in terms of 



the quantity — ^^ — • Its expression is given by: 



F(q^) = 1 + , , . 136 



In the hmit of large q, it becomes: 

%") = ^(l-6- + ...). (137) 

The dominant term at high g is a half of that one in the former calculation, Eq. (134), 
which is in relation with the fact that only the low value of k in Eq. (135) contributes to 
it, whereas in Eq. (132), the values of k around — |g*also contribute for an equal amount. 
More important here, the correction to the dominant term is of relative order q~^ instead 
of q~^. This is in agreement with an expansion in terms of powers of g~^ of the integrand 
in Eq. (135). Beyond the g"^ term which converges, the next one in q~^ diverges linearly 
and should be effectively cut-off by a factor q, hence the correction of relative order q~^. 
The same argument should apply to the full expression of the form factor, Eq. (132), but 
it turns out that the different q'^ corrections cancel out in this case, indicating that the 
way the variable q appears in this expression is of special relevance for the existence of 
g~^ corrections instead of q~^, as expected at first sight. In any case, the above discussion 
shows the difficulty to make statements in the more complicated three-body case. 

The difference as to the q^^ corrections has probably to do with the difference in 
the mathematical properties of F{(f) and F{(f) when the variable q is made complex, 
q — > |g|e*^ with 6 varying from to it. An absence of change is the sign that the function 
under consideration is a function of g^. When the above transformation is made on 
F(g^), a singularity occurs for 9 = ^n and the integration variable taking the value 

k = J'-^ k"^. This does not occur for F{q^), Eq. (133), due to the presence at the 

denominator of the term /c.g besides the quantity k'^ + ^ + K?. The weight of the singularity 
is reduced by the condition that k should be orthogonal to q. 

One can imagine that the above results be extended for some part to the three-body 
case. The absence of an exact analytic expression for the three-body wave function pre- 
vents one to make definite statements however. We nevertheless expect that corrections 
of relative order g~^ should be associated to logq factors, indicating that the results for 
the two-body system cannot be transposed as such to the three-body system. 
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